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Abstract 


We examine the Einstein equation coupled to the Klein-Gordon equation for a 
complex-valued scalar held. These two equations together are known as the Einstein- 
Klein-Gordon system. In the low-held, non-relativistic limit, the Einstein-Klein- 
Gordon system reduces to the Poisson-Schrodinger system. We describe the simplest 
solutions of these systems in spherical symmetry, the spherically symmetric static 
states, and some scaling properties they obey. We also describe some approximate 
analytic solutions for these states. 

The EKG system underlies a theory of wave dark matter, also known as scalar 
held dark matter (SFDM), boson star dark matter, and Bose-Einstein condensate 
(BEG) dark matter. We discuss a possible connection between the theory of wave 
dark matter and the baryonic Tully-Eisher relation, which is a scaling relation ob¬ 
served to hold for disk galaxies in the universe across many decades in mass. We 
show how hxing boundary conditions at the edge of the spherically symmetric static 
states implies Tully-Fisher-like relations for the states. We also catalog other “scal¬ 
ing conditions” one can impose on the static states and show that they do not lead to 
Tully-Fisher-like relations—barring one exception which is already known and which 
has nothing to do with the specihcs of wave dark matter. 



To my parents, Stewart and Carolyn. 


IV 



Contents 


Abstract [m] 

List of Tables Iviil 


List of Figures Iviiil 

Acknowledgements 

1 Introduction [T] 

1.1 General Relativity . [U 

1.1.1 Introduction to General Relativity . [T] 

1.1.2 Geometrized Units. |2] 

1.1.3 General Relativity: A Success Story. I2] 

1.2 Large-Scale Unexplained Gurvature (Dark Matter). H] 

1.2.1 Evidence for Dark Matter. H] 

1.2.2 Interpreting the Evidence . |6] 

1.2.3 Particle Theories of Dark Matter. [7] 

1.3 Wave Dark Matter. |8] 

1.3.1 Dark Matter from Geometry. IH] 

1.3.2 Wave Dark Matter as a Particle Theory . |9] 

1.4 The Baryonic Tully-Fisher Relation. |9] 

1.4.1 Observational Astronomy . [10] 

1.4.2 The Baryonic Tully-Fisher Relation. [TT] 

1.4.3 Modified Newtonian Dynamics (MOND). [TT] 

1.4.4 Dark Matter and the Baryonic Tully-Fisher Relation. [T^] 

2 The Einstein-Klein-Gordon and Poisson-Schrodinger Systems in Spher¬ 
ical Symmetry 1141 

2.1 Introduction. M 

2.2 The Spherically Symmetric Static States. [15] 

2.3 The Poisson-Schrodinger System. [TT] 

2.4 Scalings of the Poisson-Schrodinger Static States. [H] 

2.5 Approximations of the Poisson-Schrodinger Static States. [23] 


V 

























3 Wave Dark Matter and The Baryonic Tully-Fisher Relation |35] 

3.1 Introduction. [33] 

3.2 Scaling Conditions. [3l| 

3.3 Definitions. [35] 

3.4 Boundary Conditions at r = Rdm. ES 

3.4.1 Boundary Condition 1: Fixing a Length Scale. [SB] 

3.4.2 Boundary Condition 2: Fixing the Density. [37] 

3.4.3 Numerical Evidence. . m 

3.5 Central Conditions at r = 0. [37] 

3.6 Global Conditions. EH] 

3.7 Discussion and Conjectures . [10] 

3.7.1 The Slope of the Baryonic Tully-Fisher Relation. [lO] 

3.7.2 Conjectures. [H] 

A The Einstein-Klein-Gordon-Maxwell Equations [44] 

A.l Varying the Metric g . [H] 

A.2 Varying the Scalar Function /. [16] 

A.3 Varying the One-Form A . [10] 

A.4 The Einstein-Klein-Gordon-Maxwell Equations. [17] 

B The Einstein-Klein-Gordon Equations in Spherical Symmetry [IS] 

Bibliography [03 

Biography [03 


VI 


















List of Tables 


1.1 Common masses in geometrical units of time or distance. E] 

1.2 Common times in geometrical units of distance or mass. E] 

1.3 Common distances in geometrical units of mass or time. |3] 


vii 





List of Figures 


1.1 

1.2 


1.3 


Spiral galaxy NGC 4911 in the Coma Cluster. Credit: NASA, ESA, 

and the Hubble Heritage Team (STScI/AURA). 

The “bullet cluster”. In the pink regions lies the gas, which con¬ 
tains most of the baryons. In the blue regions lie the stars. Gravi¬ 
tational lensing shows that most of the mass lies in the blue regions. 
This suggests that most of the mass of clusters is non-baryonic dark 
matter. Credit: X-ray: NASA/CXC/CfA/M.Markevitch et ah; Opti¬ 
cal: NASA/STScI; Magellan/U.Arizona/D.Clowe et ah; Lensing Map: 

NASA/STScI; ESO WEI; Magellan/U.Arizona/D.Clowe et al. 

This graphic—see 28 —illustrates (a) the original Tully-Fisher rela¬ 


tion and (b) the baryonic Tully-Fisher relation. Each point represents 
a galaxy. Squares represent galaxies whose rotational velocities have 
been estimated using the width of the 21 cm line and circles represent 
galaxies whose rotational velocities have been estimated from resolved 
rotation curves. The left graph (a) plots stellar mass (which tracks 
luminosity) vs. rotational velocity, whereas the right graph (b) plots 
total baryonic mass (stellar, gas, and dust) vs. rotational velocity. 
Galaxies whose masses are gas-dominated are shown in green; they do 
not follow the original Tully-Fisher relation, but do follow the bary¬ 
onic Tully-Fisher relation. The black line has slope 4. Reproduced 
from Figure 1 in McGaugh, S. S., Schombert, J. M., Bothun, G. D., 
& de Blok, W. J. G. 2000, ApJ, 533, L99. ©AAS. Reproduced with 
permission. 
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2.2 A typical hfth excited state {n = 5) demonstrating the location of 
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(2.16) 


2.4 Five scalings of the n = 1 excited state from hg. 2.8 (shown in blue) 
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This graph shows that the approximation (2.26) is very good for large 

n and not too bad for small n . 

A ground state (n = 0) with T = 101y“^. Units are geometrized and 
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m 


eqs. (2.22) and (2.31) 


[n 


An excited state 

description. 

An excited state (n = 2). 

description. 

2.10 An excited state {n = 10). 
description. 


= 1). See the caption to hg. 2.7 for a fuller 


See the caption to hg. 2.7 for a fuller 


See the caption to hg. 2.7 for a fuller 
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1 

Introduction 


1.1 General Relativity 

1.1.1 Introduction to General Relativity 

General relativity is a physical theory written in the mathematical language of 
semi-Riemannian geometry. In this dissertation we assume a basic knowledge of 
semi-Riemannian geometry—see 31 . General relativity models the universe as a 4- 


dimensional manifold called a spacetime which is equipped with a Lorentzian metric 
of signature —+++orH -. In this dissertation we will use the —+ + + conven¬ 

tion. The fundamental equation of general relativity is the Einstein equation, which 


in geometrized units (see section 1.1.2) is 


G -1- A^f = SttT. 


( 1 . 1 ) 


Here G is the Einstein curvature tensor, a (0, 2)-tensor dehned in terms of the metric 
tensor the Ricci curvature tensor Ric, and the scalar curvature R as 


G = Ric — Rg. 

2 ^ 


( 1 , 2 ) 


The parameter A is the cosmological constant, whose effect on the evolution of the 
universe is only seen at very large scales of distance and time. In this dissertation 
we will only be studying the universe at galactic scales and take A = 0 to simplify 
matters. On the right side of eq. (1.1) we have the energy-momentum tensor T, 


another (0, 2)-tensor which encapsulates the information about the matter and energy 
content of the universe. Both G and T are symmetric tensors. 


At the highest descriptive level, eq. (1.1) says that the curvature of spacetime 
is correlated with the matter and energy content of the universe. The presence of 
matter and energy warps spacetime; conversely, the curvature of spacetime affects 
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the motion of matter and energy. Light rays (i.e., photons) and so-called “test 
particles” (objects with masses small enough that their effect on the curvature tensor 
is negligible) follow geodesics, which are purely geometrically dehned curves. Thus 
gravity is explained using geometry. 


1.1.2 Geometrized Units 


It is common in mathematical physics to choose units so that fundamental physical 
constants have magnitude 1. For example, by measuring distance in light-years and 
time in years, the speed of light is 


c = 


1 light-year 
1 year 


= iiy/yr. 


Often once this choice is made the units are omitted and it is said that “c = 1”. 

In general relativity it is customary to measure mass, time, and distance in units 
such that G = c = 1. Here G is the universal gravitational constant, not the 
Einstein tensor from eq. (1.1). Units such that G = c = 1 are called geometrized 
units. Geometrized units are useful because they simplify formulas, but they also 
lead to strange locutions. For example one can say “the Sun’s radius is 2.3 seconds” 
or “the Sun’s radius is 9.4 x 10^® kilograms” and not be talking nonsense. The hrst 
statement means that light takes 2.3 seconds to cross a distance equal to the Sun’s 
radius. An easy way to understand the second statement is to use the formula for 
the Schwarzschild radius of a black hole of mass M: 




2GM 


(1,3) 


In geometrized units, eq. (1.3) becomes 


Rsch = 2M. 


(1.4) 


Thus, an easy way to understand “the Sun’s radius is 9.4 x 10^^ kilograms” is to 
restate it as “the radius of a black hole with a mass of 9.4 x 10^^ kilograms would 
be half the Sun’s radius”. 

In tables 1^ to we have collected some common astronomical masses, times, 
and distances and list their values in geometrized units. In this dissertation we will 
adhere to the following convention: 

Convention 1. Unless otherwise noted, all formulas will be given in geometrized 
units and all physical quantities will be given in (light-)years. 


1.1.3 General Relativity: A Success Story 


The study of eq. (1.1) has yielded amazing physical and mathematical insights. In 


1916 Einstein used general relativity to explain the anomalous precession of Mercury 
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Table 1.1: Common masses in geometrical units of time or distance. 


Unit 

Seconds 

Years 

Meters 

AUs 

kilogram 

2.48 X 10-^^ 

7.85 X 10-44 

7.43 X 10-23 

4.96 X 10-39 

Astronomical Body 

Seconds 

Years 

Meters 

AUs 

Sun 

4.93 X 10^® 

1.56 X 10-43 

1480 

9.87 X 10-9 

Earth 

1.48 X lO^^i 

4.69 X 10-49 

0.00443 

2.96 X 10-44 

Moon 

1.82 X 10-^3 

5.77 X 10-24 

5.45 X 10-6 

3.65 X 10-46 

Jupiter 

4.70 X 10~9 

1.49 X 10-46 

1.41 

9.42 X 10-42 

Cygnus X-1 black hole 

7.4 X 10-^ 

2.3 X 10-42 

22 000 

1.5 X 10-^ 

Sag A* black hole 

20 

6 X 10-^ 

6 X 109 

0.04 

Milky Way 

10’^ 

10-4 

1046 

104 


Table 1.2: Common times in geometrical units of distance or mass. 


Unit 

Meters 

AUs 

Kilograms 

Solar Masses 

second 

3.00 X 103 

0.002 00 

4.04 X 1036 

203000 

day 

2.59 X 1043 

173 

3.49 X 1046 

1.75 X 1046 

year 

9.46 X 1046 

63 200 

1.27 X 1043 

6.41 X 1042 

Astronomical Time 

Meters 

AUs 

Kilograms 

Solar Masses 

age of universe 

1.31 X 1026 

8.73 X 1044 

1.76 X 1063 

8.84 X 1022 

age of solar system 

4.3 X 1026 

2.9 X 1044 

5.8 X 1062 

2.9 X 1022 


Table 1.3: Common distances in geometrical units of mass or time. 


Unit 

Kilograms 

Solar Mass 

Seconds 

Years 

meter 

1.35 X 1027 

0.000 677 

3.34 X 10-9 

1.06 X 10-46 

AU 

2.01 X 1033 

1.01 X 103 

499 

1.58 X 10-6 

light-year 

1.27 X 1043 

6.41 X 1042 

3.16 X 107 

1 

parsec 

4.16 X 1043 

2.09 X 1043 

1.03 X 10® 

3.26 

Astronomical Body 

Kilograms 

Solar Mass 

Seconds 

Years 

Sun (mean radius) 

9.37 X 1036 

471 000 

2.32 

7.36 X 10-3 

Earth (mean radius) 

8.58 X 1033 

4310 

0.0213 

6.73 X 10-46 

Moon (mean radius) 

2.34 X 1033 

1180 

0.005 79 

1.84 X 10-46 

Jupiter (mean radius) 

9.41 X 1034 

47300 

0.233 

7.39 X 10-9 
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and predicted the correct angle at which light wonld be deflected by the snn. He 
also predicted the gravitational redshift of light, which was definitively observed in 
1959, the existence of gravitational waves, for which indirect evidence was discovered 
in 1993, and the possibility of gravitational lensing, which is cnrrently being nsed 
extensively to map the distribntion of matter in the nniverse. The Schwarzschild 
solntion, discovered in 1916, snggested the possibility of black holes, whose existence 
is now accepted. 

In addition to these snccesses in solar system dynamics and astrophysics, general 
relativity is indispensible in cosmology. The Friedmann-Lemaitre-Robertson-Walker 
metric provides the standard model of the nniverse and includes the Big Bang. Also, 
the most popular way to explain the accelerating expansion of the universe observed 
in 1998 is with the cosmological constant A. This often goes under the name of “dark 
energy”. 

Fundamental theoretical work in geometry contributing to our understanding of 


the universe includes the positive mass theorem first proved by Schoen and Yau 39 
40 and later Witten 51 , the singularity theorems of Hawking and Penrose 33, 16 


and the proofs of the Riemannian Penrose inequality by Huisken and Ilmanen 
and Bray [5]. 


18 


1.2 Large-Scale Unexplained Curvature (Dark Matter) 

1.2.1 Evidence for Dark Matter 

In 1933, the astronomer Fritz Zwicky measured the velocites (via gravitational red- 


shift) of galaxies belonging to the Coma cluster (see fig. 1.1). The galaxies seemed to 
be moving too fast; that is, given the total mass of the cluster implied by the amount 
of visible matter, the gravitational pull on most galaxies should have been nowhere 
near strong enough to keep the galaxies from flying off into deep space, given their 
velocities. In fact, a quantitative analysis showed that for the galaxies to remain 
gravitationally bound, the total mass of the Coma cluster must be at least ten times 
the mass of the visible matter. Zwicky hypothesized that “Dunkle Materie”—dark 
matter—was present in the Coma cluster, but invisible because it did not interact 
with light- 


-see 


52 


In the decades since Zwicky’s proposal, other evidence has arisen that supports 
the hypothesis of dark matter. Constraints from Big-Bang nucleosynthesis point 
to a baryon density only 5% of the critical density, whereas detailed studies of the 
anisotropies in the cosmic microwave background (CMB) indicate that the universe 
is flat (i.e. at critical density). Furthermore the anisotropies in the CMB are so small 
that without dark matter to “seed” large-scale structure growth, the universe could 
not be as far along in its evolution as it is today. These arguments suggest that most 
of the matter/energy content of the universe is non-baryonic. 

Gravitational lensing studies have also produced evidence for dark matter. Clus¬ 
ters of galaxies can act as giant gravitational lenses and produce multiple images of 
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Figure 1.1: Spiral galaxy NGC 4911 in the Coma Cluster. Credit: NASA, ESA, 
and the Hubble Heritage Team (STScI/AURA). 


background galaxies. This is an example of “strong gravitational lensing”. Strong 
gravitional lensing can provide an estimate of the mass of the lens, and the evidence 
consistently points to the mass of clusters being much greater than the visible mass. 
Studying “weak gravitation lensing” is also very popular. One famous study 24 of 
the “bullet cluster” describes observations of the collision of two clusters of galaxies 
(see £g. 1.2). In clusters, most of the baryons do not he in stars but in the intra¬ 
cluster medium (ICM) in the form of hot gas. In a collision of clusters, it is only the 
gas that really collides and is slowed by friction. The stars in the galaxies are much 
too far apart to collide and pass through each other relatively unaffected. Thus, a 
collision between clusters separates the ICM from the stars. Figure 1.2 shows this 
separation. In pink is the ICM, and on the right is a very obvious bow shock from 
which the epithet “bullet cluster” derives. In blue are the galaxies with their stars. 
Remarkably, gravitational lensing reveals that the blue regions contain much more 
matter than the pink regions, exactly the opposite of what we would expect if all the 
matter were baryonic. The bullet cluster is usually cited as one of the best pieces of 
evidence for the existence of large quantites of non-baryonic dark matter which can 
be separated from baryonic matter. 

The evidence for dark matter that is most relevant for this dissertation comes from 
studying disk galaxies. In the 1970s, Vera Rubin and Kent Ford began publishing 
“rotation curves” for disk galaxies which plotted the rotational velocity of stars, 
gas, and dust in a galaxy as a function of radius- 


-see 


36 35 . If all matter were 


baryonic, rotation curves would decline signihcantly with increasing radius, and stars 
at the outer edges of a galaxy would orbit much more slowly than stars close to the 
center. Instead, most rotation curves are roughly flat, indicating that most baryonic 
matter in a disk galaxy orbits at the same speed. The standard explanation for this 
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Figure 1.2: The “bullet cluster”. In the pink regions lies the gas, which 
contains most of the baryons. In the blue regions lie the stars. Gravita¬ 
tional lensing shows that most of the mass lies in the blue regions. This sug¬ 
gests that most of the mass of clusters is non-baryonic dark matter. Credit: 
X-ray: NASA/CXC/CfA/M.Markevitch et ah; Optical: NASA/STScI; Magel¬ 
lan/U.Arizona/D. Clowe et ah; Lensing Map: NASA/STScI; ESO WFI; Magel- 
lan/U.Arizona/D.Clowe et ah 


phenomenon (but see section 1.4.3 below) is that each disk galaxy is embedded in 
a large spherical “halo” of dark matter which dominates the baryonic matter by a 
factor of 5 or more and which flattens out the rotation curve. 


1.2.2 Interpreting the Evidence 

We have described many of the types of observations that lead most experts to believe 
in the existence of dark matter. It is important to understand that dark matter has 
never been directly observed. All astronomical observations ultimately originate 
with the act of gathering electromagnetic radiation with telescopes. Because dark 
matter is, almost by dehnition, matter which does not interact with light, it could 
never be observed astronomically, only inferred. For this reason, arguments about 
dark matter which rely mainly on astronomical observations are especially “theory¬ 
laden”. For example, at the most basic level all the arguments from the previous 
section assume that the theory of general relativity and its low-held, non-relativistic 
limit (Newtonian gravity) are basically correct. If this assumption is not justihed, 
then the reason our observations don’t make sense is that we just don’t understand 
gravity, and the arguments for dark matter fall hat. Some physicists have pursued 
this line of thought and proposed various modihcations of general relativity and 
Newtonian gravity. The most popular is Modihed Newtonian Dynamics, or MOND. 


We briehy discuss MOND in section 1.4.3 
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Even if general relativity and Newtonian gravity are basically correct, there is still 
a problem of interpretation of the data. In general relativity, gravity is a manifesta¬ 
tion of the curvature of spacetime. The mysteries which the dark matter hypothesis 
is supposed to solve are gravitational mysteries, or, we could just as easily say, cur¬ 
vature mysteries: we observe curvature of spacetime on large scales for which we 
are unable to account. Whether this anomalous curvature is due to some kind of 
particle which we could detect with sufficiently advanced technology or is a purely 
geometrical phenomenon is an open question. 

As an analogy, consider the cosmological constant term Ag which appears in 


Einstein’s equation (1.1). We have the freedom to place it on the left or right side 
of the equals sign without changing the mathematics. On the one hand if it is 
on the right, then it contributes to the energy-momentum tensor T which contains 
the information about the matter and energy content of the universe. Thus we are 
tempted talk of “dark energy” which permeates spacetime. On the other hand if it 
is on the left, then we can incorporate it in the Einstein curvature tensor and we are 
tempted to talk of a cosmological constant, a purely geometrical phenomenon which 
affects how matter and energy curve spacetime. 

Similarly, the large-scale unexplained curvature we observe might be due to large 
quantities of dark matter particles which are curving spacetime. Or, the curvature 
might be more embedded in the laws of physics than physical stuff embedded in the 
universe. 


1.2.3 Particle Theories of Dark Matter 

The most popular theories of dark matter today postulate a dark matter particle in 
the ACDM cosmological paradigm. A is the cosmological constant and the acronym 
CDM stands for “cold dark matter”. Cold dark matter is dark matter which travels 
at non-relativistic speeds, as opposed to hot dark matter which travels at relativistic 
speeds. It is generally agreed that hot dark matter is not viable because it does not 
lead to the type of structure formation necessary to create our universe. Cold dark 
matter is by far the most popular type of dark matter under consideration. The 
most popular particle candidates which are cold are WIMPs. 

Weakly Interacting Massive Particles (WIMPs) interact with themselves and 
other particles via the gravitational force and the weak force. The inspiration for 
WIMPs comes from particle physics. Some simple extensions of the Standard Model 
involving supersymmetry predict WIMP-like particles, and it is a natural idea to 
identify these hypothetical particles with the mysterious dark matter. However, the 
theory of WIMPs does have its problems. The odds against supersymmetry are 
increasing as the Large Hadron Collider continues to operate without discovering 
evidence for it, so the main theoretical structure supporting WIMPs is not looking 
strong. Also, numerical simulations of WIMPs in galaxies continue to support the 
notion that there should be “cusps” of WIMPs at the centers of galaxies—that is, 
most every numerical simulation produces a distribution of WIMPs in which the 
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density at the galactic center is nnbonnded. Astronomical observations which at¬ 
tempt to map the distribntion of dark matter in actnal galaxies do not observe cnsps 
bnt “cores”—the density of dark matter levels off at the center. This is known as 
the cnsp/core problem. Finally, direct detection experiments for WIMPs have been 
operating on Earth for decades and have never reported any nnambignonsly positive 
resnlts. 

There are many, many hypothetical particles which theorists have pnt forward as 
candidates for dark matter. For a good introdnction to the possibilities we recom¬ 
mend 


21, 30 


1.3 Wave Dark Matter 


1.3.1 Dark Matter from Geometry 

The theories of dark matter described in the previons section draw on particle physics 
for their inspiration. Bnt since the dark matter mystery is fnndamentally a gravi¬ 
tational mystery, and given the nnmerons snccesses of general relativity as a gravi¬ 


tational theory (section 1.1.3), perhaps we should turn to general relativity for in¬ 
spiration. Moreover, if we take seriously the key insight of general relativity that 
gravity is geometry, then perhaps we should turn to geometry for inspiration. In 
2013 Hubert Bray published a paper in which he proposed to examine what hap¬ 
pens if we do not make the usual assumption that the connection of the universe is 
the Levi-Civita connection. He found that the deviation of a spacetime connection 
from the Levi-Civita connection could be described using a scalar function and that 
this function would obey the Klein-Gordon wave equation. The full equations of the 
theory are the coupled Einstein-Klein-Gordon equations 



df ® df + df ® df 

'f2 




(1.5a) 

(1.5b) 


The reader will note that eq. (1.5a) is just the Einstein equation (1.1) with a par¬ 


ticular form for the energy-momentum tensor. Equation (1.5b) is the Klein-Gordon 


equation and / is the scalar function which describes the deviation of the connec¬ 
tion from the Levi-Givita connection. The parameter T is a fundamental constant 
of the theory. Ghapter of this dissertation disusses some results concerning the 
Einstein-Klein-Gordon system in spherical symmetry. 

Bray proposed that the scalar function / should be identihed with dark matter 
and sought to model dark matter in spiral galaxies. He obtained intriguing results 
with respect to the growth of spiral structure. In a subsequent paper he ob¬ 
tained results indicating that the wave dark matter is consistent with the existence 
of “shells” in elliptical galaxies. In 32 , Alan Parry examined wave dark mat¬ 
ter in the context of dwarf spheroidal galaxies. In this dissertation we propose to 












investigate wave dark matter in the context of disk galaxies and its prospects for ac¬ 
counting for an astronomical phenomenon called the baryonic Tully-Fisher relation. 


We introduce the baryonic Tully-Fisher relation in section 1.4 


Equations (1.5a) and (1.5b) (along with Maxwell’s equations for electromag¬ 


netism) can be derived using a variational principle from the following action: 


H9J.A) = 


\ 


R-2A- Wtt 


Wl 

Y‘2 


+ 1 / 1 ^+ 


dV. 


( 1 . 6 ) 


Here H is a one-form, the electromagnetic potential. In appendix we use a varia¬ 
tional principle to derive the Einstein-Klein-Gordon-Maxwell equations. The action 


(1.6) is remarkably simple, and, if the theory of wave dark matter is correct, describes 


most of the matter and energy content of the universe for most of its history. For in 
its earliest times the universe was dominated by electromagnetic radiation, at later 
times by dark matter, and today is dominated by dark energy. 


1.3.2 Wave Dark Matter as a Particle Theory 

Although Bray motivated the theory of wave dark matter using a geometrical per¬ 
spective from general relativity, the theory has been investigated before under other 


names such as scalar held dark matter (Isj pGj [4^ |T^ S 

iO and boson star dark 

matter or Bose-Einstein condensate (BEG) dark matter 

43 

19 

22 

25 

42 

48 . 


Results concerning axionic dark matter are also relevant as the axion is also a scalar 
held particle. Under all these other names the motivation is a more conventional one 
from particle physics. But the underlying equations—the Einstein-Klein-Gordon 
equations (1.5a) and (1.5b)—are the same, so which motivation we prefer is irrel¬ 


evant for the theory’s predictions. In a particle physics approach, the relationship 
between the fundamental constant T and the particle mass m is 


m = — = 2.09 X 10" 
c 


23 


eV 


T 


lly-i 


(1.7) 


If wave dark matter turns out to be the correct theory of dark matter, the question 
of which approach corresponds more closely to reality might remain ambiguous. The 
only way we foresee this question being resolved is the direct detection of the dark 
matter particle. 


1.4 The Baryonic Tully-Fisher Relation 

Ghapter of this dissertation discusses the relevance of wave dark matter for an 
empirical astrophysical relation concerning disk galaxies known as the baryonic Tully- 
Fisher relation. In this section we introduce this relation. 
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1.4-1 Observational Astronomy 


We first need to introduce a few basic concepts from observational astronomy. 

The luminosity of a galaxy is a measure of the amount of energy emitted by the 
galaxy per unit time. It is correlated with brightness, but brightness usually refers 
to the galaxy’s output in the visible wavelengths, whereas luminosity accounts for 
the galaxy’s output at all wavelengths. It is a distance-independent quantity. The 
baryonic mass of a galaxy is the total mass of the galaxy’s baryons, i.e. the ordinary 
matter of the periodic table. Measuring/estimating luminosity and baryonic mass 
is a very complicated business which we will not attempt to describe. Note that 
both concepts suffer from being rather ill-dehned because galaxies do not have sharp 
edges. Different conventions for how to dehne where a galaxy ends can lead to wildly 
different values for these parameters in the literature. 

Emission lines are sharp peaks in the galaxy’s electromagnetic spectrum which 
correspond to various transitions between electron energy levels in atoms. Their 
presence or absence can reveal a lot about the galaxy’s composition. For example, 
in neutral hydrogen atoms there is a “hyperhne” transition which occurs very rarely. 
When it occurs, a photon is emitted with a wavelength of ~21 cm at a frequency of 
~1420MHz. A source whose spectrum has a 21cm line contains large amounts of 
neutral hydrogen (“large amounts” because the transition is so rare). The spectra 
of galaxies always contain the neutral hydrogen line because galaxies contain huge 
“H I” regions of neutral hydrogen gas. This line is particularly useful in observational 
astronomy because the microwave wavelength penetrates cosmic dust whereas visible 
wavelengths do not. 

Emission lines are subject to the Doppler effect. If a source is moving towards 
us its spectrum is shifted towards shorter wavelengths/higher frequencies and if it 
is moving away from us its spectrum is shifted towards longer wavelengths/lower 
frequencies. These shifts are most commonly called “blueshift” and “redshift”, re¬ 
spectively, because in the visual spectrum the color blue lies on the low wavelength 
end and the color red on the high wavelength end. A disk galaxy’s spectrum is the 
composite of the spectra of its stars, gas, and dust. If we are observing it at an 
angle, then because it is rotating the components will have slightly different Doppler 
shifts. This causes Doppler line broadening of the emission lines. The width of the 
emission lines then provides some measure of the average rotational velocity of the 
galaxy. Galaxies which rotate faster have broader emission lines. 

The rotation curve of a disk galaxy gives the rotational velocity of the galaxy as 
a function of radius. We have already introduced this concept in the discussion of 
the evidence for dark matter in section 1.2. 1[ Most disk galaxies have rotation curves 
which are roughly flat. Obtaining a galaxy’s rotation curve involves analyzing the 
Doppler shifts in the galaxy’s spectrum at different radii. 
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1.4-2 The Baryonic Tully-Fisher Relation 


In a seminal paper 47 , Richard Tully and James Fisher reported that for disk 
galaxies luminosity L seemed to be related to the width of the 21cm line w by a 
power law: 

L a w^. 

Because the width of the 21 cm line is a proxy for a galaxy’s rotational velocity, 
qualitatively the Tully-Fisher relation says that brighter galaxies rotate faster. This 
makes sense because a galaxy’s brightness is correlated with its mass—more mass 
means more stars and more light, and also stronger gravity and faster rotation. 

McGaugh proposed in an update to the Tully-Fisher relation. He noted 
that many smaller galaxies (low surface brightness galaxies, or LSBs) fell below the 
Tully-Fisher line in luminosity vs. line width log-log space because much of their 
mass was in the form of dim gas and dust instead of luminous stars. The Tully-Fisher 
relation, while linear for large galaxies, broke down for these LSBs. See £g. 1.3 The 


problem was that the luminosity measure was effectively counting only the stellar 
mass and ignoring the mass of gas and dust—and rotational velocity depends on 
mass irrespective of its form. McGaugh noted that when the mass of gas and dust 
was included and total baryonic mass was plotted vs. rotational velocity, a linear 


relationship was restored over the entire range of galactic masses. Again, see £g. 1.3 
The baryonic Tully-Fisher relation has the form 


Mb a 


( 1 , 8 ) 


where Mb is baryonic mass and v is rotational velocity. For galaxies whose rotation 
curves can be resolved, v can be the maximum observed velocity, or the average 
velocity, or the velocity at a hxed radius. Since the rotation curves are close to flat 
all of these choices are effectively equivalent. In studies where rotation curves are 
not available v is obtained via measuring the neutral hydrogen line width. 

The slope x of the baryonic Tully-Fisher relation has been the subject of many 
investigations (see [^)- The literature seems to agree 

that 3 < a: < 4 with numbers closer to 4 being favored, but beyond that there is 
little agreement. Most of the uncertainty seems to stem from disagreement about 
how to estimate Mb- In addition, there is a danger of bias because a slope of 3 is 
thought to be more consonant with the standard AGDM cosmology whereas a slope 
of 4 is consonant with Modified Newtonian Dynamics, the most popular alternative 


to the dark matter hypothesis (see section 1.4.3). 


1-4.3 Modified Newtonian Dynamics (MOND) 


Recall the discussion in section 1.2.2 in which we pointed out that in trying to 


account for the confusing observations described in section 1.2.1, we can either seek 
a solution within the paradigm of general relativity and Newtonian gravity or look 
for a new theory of gravity. The most popular proposed solution of the latter type is 
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known as Modified Newtonian Dynamics or MOND for short 29, 8j. MOND, while 


it has other issues, can claim to explain the flat rotation curves of spiral galaxies 
and the baryonic Tully-Fisher relation. Indeed, it was designed for this purpose. In 
essence, whereas the combination of Newton’s second law and law of gravity gives 
an acceleration due to gravity 

GM 


a = 


MOND postulates 


a = 


GMoiq 


for an acceleration a much less than a threshold acceleration oq. The inclusion of the 
threshold acceleration is to leave solar system dynamics virtually unchanged. One 
immediately sees that for circular motion where a = v^/r, we get v = {GMaoY^^ 
(velocity independent of radius) and M/v'^ = {Gao)~^ (a Tully-Fisher relation), which 
seems promising. However, the theory has its own conflicts with data. One of the 
most problematic is that although MOND was created to get rid of the missing 
mass problem in galaxies, it has a missing mass problem at the level of clusters 10 
38 . Even more problematic is the “bullet cluster” described in section 1.2.1, whose 


existence seems to demonstrate that dark matter exists in large quantities and can 
be separated from baryonic matter. MOND remains a minority viewpoint among 
astrophysicists. Even so, there is still the important question of why it works so well 
for disk galaxies. 


1.4-4 Dark Matter and the Baryonic Tully-Fisher Relation 

We have introduced the baryonic Tully-Fisher relation because it may have relevance 
for discriminating between theories of dark matter. Dark matter comprises the bulk 
of a galaxy’s mass, is responsible for the flatness of the rotation curves, and is the 
dominant factor determining the v in the relation (1.8). However the mass that 


appears in the relation is the baryonic mass, not the total mass or the dark mass. 
Thus the baryonic Tully-Fisher relation gives a quantitative link between the dark 
and the baryonic mass. A successful theory of dark matter needs to be consistent 
with the existence of the relation. We will have more to say on this subject in 
chapter after we have described our results. 
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Figure 1.3: This graphic—see 28 —illustrates (a) the original Tully-Fisher relation 


and (b) the baryonic Tully-Fisher relation. Each point represents a galaxy. Squares 
represent galaxies whose rotational velocities have been estimated using the width 
of the 21 cm line and circles represent galaxies whose rotational velocities have been 
estimated from resolved rotation curves. The left graph (a) plots stellar mass (which 
tracks luminosity) vs. rotational velocity, whereas the right graph (b) plots total 
baryonic mass (stellar, gas, and dust) vs. rotational velocity. Galaxies whose masses 
are gas-dominated are shown in green; they do not follow the original Tully-Fisher 
relation, but do follow the baryonic Tully-Fisher relation. The black line has slope 4. 
Reproduced from Figure 1 in McGaugh, S. S., Schombert, J. M., Bothun, G. D., & 
de Blok, W. J. G. 2000, ApJ, 533, L99. (c)AAS. Reproduced with permission. 
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2 


The Einstein-Klein-Gordon and 
Poisson-Schrodinger Systems in Spherical 

Symmetry 


In this chapter we study the Einstein-Klein-Gordon system in spherical symmetry 
from a mathematical perspective. In the next chapter we will consider the relevance 
of our results for theoretical astrophysics and the problem of dark matter. 


2.1 Introduction 


The Einstein-Klein-Gordon system couples the Einstein equation (1.1) together with 
the Klein-Gordon equation for a complex-valued function. The full system is 


G -I- A^f = Svr 

□/ = TV- 


df®df + df®df ( \dfy 


^-2 


^-2 


+ \fn9 


(2.1a) 

(2.1b) 


Here T is a positive constant. This system can be obtained via a variational principle 
from the action 


= PR- 2A - levr 


m 

^-2 


+ l/n dV- 


Requiring the metric g and the scalar function / to be critical points of the action 
leads to eqs. (2.1a) and (2.1b). See appendix [A} 

A full solution to the system (113 consists of a spacetime N with a Lorentzian 
metric g and a scalar function function f : N C such that eqs. (2.1a) and (2.1b) 


14 










are satisfied. We are interested in studying spherically symmetric static solutions to 
the Einstein-Klein-Gordon system where N has the topology of As an ansatz for 
the metric g we take 

1 _ 2M{r) \ ^^2 ^2.2) 

r J 


9 


= dt^ + 


Here M (r) and V(r) are functions of the coordinate r only and dO"^ + r‘^ sin^ 6 dcj)'^ is 
the standard area form on the coordinate sphere of radius r. The reason for choosing 
this particular ansatz, with the functions M(r) and V{r), is that these two functions 
have natural physical interpretations—in the Newtonian limit, M(r) will turn out to 
be the total mass contained inside the sphere of radius r and V (r) will turn out to 
be the Newtonian gravitational potential. For clarity we dehne 


<h(r) = 1 - 


2M(r) 


so that the metric (2.2) can be written as 


g = — dt^ + <h(r) ^ dr^ + dO^ + sin^ 6 dcj)^. 


(2.3) 


(2.4) 


2.2 The Spherically Symmetric Static States 


The scalar function / must be spherically symmetric, a function of t and r only. We 
make the additional ansatz that it can be written in the form 


f{t,r) = F{r)e^-^ 


(2.5) 


for some real function F{r) and real number u. A solution in which / has the form 
) is called a static state. In appendix]^ we show that with these ansatze, solving 
the Einstein-Klein-Gordon system reduces to solving the following system of coupled 
ordinary differential equations: 



Mr = ^ [(T^ + F^ + <!>F^] (2.6a) 

^ ^ F^ - 4)^2] (2.6b) 

Frr + -Fr + WF, + (T^ - F. (2.6c) 

Note that the dependence on t has disappeared, which is why solutions to these 
ODEs are called static states. Our notation for a solution will be 


(a;;M,K,F) 


because to specify a solution we need to specify the frequency oj appearing in eq. (2.5) 
and the three functions M, V, and F. Solutions can be found by numerical integra¬ 
tion using a computer. 
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For initial conditions we take 


M(0) = 0 (2.7) 

V^(O) = K) (2.8) 

F(0) = Fo > 0 (2.9) 

F^(0) = 0. (2.10) 


Equation (2.7) must hold for smoothness of the solution at r = 0; the function 
<F(r) cannot become inhnite as r 
eq. 


( 2 . 8 ) 


0. We do not place any restriction on Vq in 
-but see Convention below. In eq. (2.9|) we require that Fq be positive; 


this is because if (w; M, V, F) is a solution, so is (ca; M, V, —F), and thus we might as 
well take F(0) ^ 0. We exclude F(0) = 0 because this leads to the trivial solution. 


Finally, eq. (2.10) is another requirement for smoothness at r = 0 because of the 


spherical symmetry. 

Note that we have the freedom to add an arbitrary constant V to the potential 


function V{r). Looking at the form of the ODEs (2.6), we see that if (w; M, V, F) is 
a solution, then so is M, V + V, F). This corresponds to just a rescaling of the 
t coordinate by a factor of in the spacetime metric (2.4). Thus, adding a constant 


to V(r) amounts to a change of coordinates which does not affect the solution. 
For brevity we make the following dehnitions: 


Definition 1. 


Moo = lini M(r) (2.11) 

r—>00 

Ko = lim V(r). (2.12) 

r^oo 


Numerical experimentation with values for cj, Vo, Fq reveals that almost all solu¬ 
tions blow up as r —>■ oo; that is. Moo = Voo = +oo and limr^oo F(r) = +oo. We are 
not interested in these solutions that blow up because they cannot represent any¬ 
thing physical. We are interested in the rare solutions where Moo and Voo are hnite 
and limr .^00 F(r) = 0. In chapter we will use these special solutions to model dark 
matter in disk galaxies. In a moment we will begin to describe the characteristics of 
these solutions, but hrst we introduce the following convention: 

Convention 2. For a solution (cu; M, V, F) with Voo hnite, we will assume that 


Ko = 0. 


(2.13) 


Any solution with Voo finite can always be made to satisfy eq. (2.13) by per¬ 
forming the transformation mentioned above regarding the function E(r); this will 
involve changing the value of u, but as we said above, this is really just a change of 
coordinates. There are at least two good reasons for using this convention: (1) The 
metric (2.4) is then asymptotic to Minkowski spacetime at infinity. (2) Since the 
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physical interpretation of V(r) is as the Newtonian potential, this convention agrees 
with the common convention in physics to have the potential equal zero at inhnity. 

With Convention!^ in place we now describe the special solutions just mentioned. 
(From now on when we will use “static state” to refer to one of these special solu¬ 
tions.) The static states come in the form of ground states and excited states. In 


£g. 2.1 we have graphed a ground state and hrst three excited states with T = 100 


and uj = 99.9. For fixed T and u, there are countably many solutions corresponding 
to n = 0,1, 2 ... where n is the number of zeros (or “nodes”) of F. The static states 
live in a 3-parameter space, where two of the parameters are continuous and one is 
discrete (for example, T, cj, and n —but there are other ways of parametrizing the 
static states.) 


2.3 The Poisson-Schrodinger System 


To gain more insight into the static states, we consider the system (2.6) and make 
the following approximations: 


2V 

e sa 


$ 


K 


T\\Vl 


0 , 




T 


0 


^ 1 
Y ^ h 


T\\Fl 


0 . 


(2.14) 

(2.15) 


The approximations in (2.14) may be interpreted as saying that the metric (2.4) 


is close to the Minkowski metric (the low-held limit), and the approximations in 


(2.15) may be interpreted as saying that the group velocities of wave dark matter 


are much less than the speed of light (the non-relativistic limit). Applying all these 


approximations to the system (2.6) leads to the system 

Mr = • 2F^ 




(Frr + ~Fr 

2T V r 


= (T - a; + TV)F. 


(2.16a) 

(2.16b) 

(2.16c) 


This system is the Poisson-Schrodinger system written in spherical symmetry. Here 
we have a special case of the well-known fact that the Poisson-Schrodinger system 


is the non-relativistic limit of the Einstein-Klein-Gordon system—see 11 . Also see 


kg. 2.3 Thus, to understand the system (2.6) in the low-held, non-relativistic limit 


it suffices to understand the system (2.16) 


Equations (2.16a) and (2.16b) immediately lead to the physical interpretations 


of M{r) as the total mass contained inside a sphere of radius r and V{r) as the New¬ 
tonian gravitational potential. The quantity 2F‘^ represents the mass-energy density 


p of the static state. Combining eqs. (2.16a) and (2.16b) gives the Poisson equation 
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Figure 2.1: The ground state (top row) and first three excited states (second, third, 
and fourth rows) with T = 100 and u = 99.99. The vertical black line in each plot 
marks the location of i^oM- 
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AV = Airp in spherical symmetry. For in spherical coordinates, the Laplacian of V 
is given by 14-r + (2/r)14-; thus 


M 


2 M M 


A.l/ = — +- =^=A7v 2F^ = Anp. 




The third equation (2.16c) is the most interesting. It is the Schrodinger half of the 
Poisson-Schrodinger system. On the left side we recognize the Laplacian of F. Let 
us consider for a moment the differential equation 


Frr H- Fr = kF. 

r 


(2.17) 


If we let h{r) = rF{r), then the corresponding differential equation for h is hrr = kh. 
Thus the general solution to eq. (2.17) is 


F{r) = \ A + - 


. ^Vkr — 

-h B- - 


k>0 
k = 0 

^^ sin(-F^r) cos(V^r) ^ ^ Q 


(2.18) 


Thus, on any interval where the expression k{r) = 2T(T—a;-l-Tl/(r)) is roughly con¬ 


stant, the solution to eq. (2.16c) will look like one of the three solutions in eq. (2.18) 


More explicitly, when k{r) is negative, F will exhibit oscillatory behavior, and when 
k{r) is positive, F will exhibit exponential behavior. Since V is an increasing func¬ 
tion, k{r) increases with r. Remembering that one of our requirements on a static 
state solution {u; M,V, F) is that Moo < oo, looking at eq. (2.16a) we see that F 


must exponentially decay after a certain point. Since we are also requiring ITo = 0, 
we must take a; < T so that limj-^oo = 2T(T — u) >0. On the other hand, 
k{r) cannot be positive (i.e., F cannot have exponential behavior) for all r; looking 


at eq. (2.18), we see that this would be in contradiction with the requirements from 


spherical symmetry that Mr{0) = Fr{0) = 0. Thus, the only situation consistent 
with our requirements is when 0 < ca < T, such that k{r) begins negative and lim¬ 
its to 2T(T — ca). The corresponding behavior of F is to start out oscillating and 
then to switch over to exponential behavior. The switch occurs at the point where 


k{r) = 0; we label this point Rbm (see £g. 2.2) in anticipation of the role it will play 
in chapter as the effective radius of a dark matter halo. 

Definition 2. Given a static state {oj] M,V, F), we dehne Rbm to be the radius at 


which the function F switches from oscillatory to exponential behavior. See £g. 2.2 


The previous discussion concerned the Poisson-Schrodinger system, but a similar 
discussion applies for the Einstein-Klein-Gordon system, in which case the expression 
which controls the oscillatory/exponential behavior of F is k{r) = T‘^. Thus 
for the Einstein-Klein-Gordon system, Rbm satisfies — cj2g-2y(i(;DM) = which 
implies 

^g-n«DM) = Y. (EKG) (2.19) 
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Figure 2.2: A typical fifth excited state {n = 5) demonstrating the location of 
Rdm (see dehnition 1^. We have omitted the plot of the potential V{r). To the 
left of Rdm, F{r) exhibits oscillatory behavior. To the right of Rbm, F{r) exhibits 
exponentially decaying behavior. 


For the Poisson-Schrodinger system, Rbm satishes T — oj + V (i?DM) = 0, which 
implies 

F(i?DM) = o; - T. (PS) (2.20) 

Solving for one of the Einstein-Klein-Gordon spherically symmetric static states 
on the computer is conceptually easy but computationally intensive. If we hx T and 
cj, then to solve for the static state of order n we proceed as follows: Fix a value 
for Vq consistent with k{0) < 0 and then begin varying Fq. Increasing Fq means 
that we get more oscillations before reaching Rbm and decreasing Fq means fewer 
oscillations. Thus to get a ground state or a particular excited state we can only 
consider values for Fq in a particular range. In that range, there is only one value of 
Fq such that F(r) exponentially decays for all r > Rbm- For all other values of Fq, 


the solution includes the exponential growth term from eq. (2.18) which dominates 
as r —> 00 , and we have agreed to ignore these solutions. Once we have found (to 
within the desired accuracy) the correct value of Fq, we then consider the fact that 
we probably do not have Ko = 0. We then vary Vq and begin again. By varying Vq 
in an outer loop and Fq in an inner loop, we can End a solution (a;; M, V, F) of order 
n. 

Next we consider the effect of changing uj. Recall that physically, T is to be 
regarded as a fundamental constant of nature. On the other hand, a; is a parameter 
we can vary. Once we have chosen u, there is then a unique ground state, first excited 


state, second excited state, etc. We find experimentally using the exact system (2.6) 


that for values of u close to T we are in the low-field, non-relativistic limit, but as u 
decreases, the mass of the state increases, Rdm decreases, and we leave the low-field. 
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non-relativistic limit. This is illustrated in fig. 2.3 for the ground state. We repeat 
again that the static states live in a 3-parameter space, where two of the parameters 
(e.g. T,a;) are continuous and one (e.g. n) is discrete. 

In the next two sections we discuss properties of the Poisson-Schrodinger static 
states. The reader should keep in mind that all these properties are valid for the 
Einstein-Klein-Gordon static states provided they are in the low-field, non-relativistic 
regime. 


2.4 Scalings of the Poisson-Schrodinger Static States 


For the Poisson-Schrodinger system, it is a remarkable fact that up to scalings, there 
is a unique state for each n ^ 0. This statement is made precise in the following 
theorem. 


Theorem 1. Suppose (ca; M, V, F) solves the Poisson-Schrodinger system (2.16) for 
a particular value of T. Let a,/3 > 0. Then {dj;M,V,F) defined by 


r = a ^(3 

(2.21a) 

M = ap-^M 

(2.21b) 

V = a^/3-^V 

(2.21c) 

F = a^F 

(2.21d) 

T = 

(2.21e) 

(T-(h) =a^{T -oj) 

(2.21f) 


is also 
Proof. 


a solution of (2.16) with T replaced by T. 
With the above definitions. 


M(f) = a(3~^M{a(3f) 
V{f) = a^ft~‘^V{al3f) 
F{f) = a^F^ajdf). 


Therefore 


Mf = a(d~^Mr{a/3f) ■ aft 
= a‘^ft~‘^Mr{aftf) 

= a^ft~‘^ ■ Sn^aftrYF{aftfY 
= S'Kf‘^{a^F{a(3fY) 

= 87rf‘^F^ 
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g; = 99 





g; = 85 _^^_ 

Figure 2.3: Solving for the ground state with T = 100 and u = 99,95,90,85. As 

u decreases we leave the low-field, non-relativistic limit. The more solid graphs are 


the solutions obtained using the exact Einstein-Klein-Gordon system (2.6), and the 


fainter graphs are the solutions obtained using the low-field, non-relativistic limit 


Poisson-Schrodinger system (2.16). 
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and we see that M{r) satisfies eq. (2.16a). Similarly 


Vf = 0^(3 ‘^Vr{a/3r) ■ a/3 
= a^/3~^Vr{a/3f) 


= a'^f3 


3 1 M{a/3r) 


{aj3ry 

a{3~^M{a{3f) 


f2 


M 

J.2 


and we see that V{r) satisfies eq. (2.16b). Finally 

1 


2(3^ 


a^Frr{aj3r) ■ iaj3Y + -a‘^FJaj3r) ■ ia(3) 

r 


2T 


= a‘^ ■ — i Frria/3r) H- —Fr{a(3r) 


a/3r 


= a\T -u + TV{a/3f))F{a(3f) 

= {a\T -u)+ (3^Ta^/3-^V{a/3f)) ■ a^F{a/3f) 
= {t-0 + tV)F 


and we see that F{r) satisfies eq. (2.16c). 


□ 


Note that (3 effectively controls the fnndamental constant of natnre T and a can 
be nsed to control the “peakedness” of the solntion. States which are more compact 
have higher mass; states which are more spread ont have lower mass. In fig. |2.4| we 


show the effect of scaling a gronnd state nsing valnes of a near 1. It follows from 
these scalings that having fixed T, there is a one-parameter family of static states of 
order n. Once we fix a farther characteristic of the solntion, a static state is nniqnely 
determined for each n, giving ns a sequence of static states. We call this “fixing a 
scaling” or “imposing a scaling condition”. 


2.5 Approximations of the Poisson-Schrodinger Static States 

We can approximate a Poisson-Schrodinger static state around Rbm by using a linear 
approximation for V{r): 


V{r) ^ P(itDM) + V\RnM)ir - Rdm) 


— V (Rbm) + 


M{R 


dmJ 


C)2 

■^DM 


(r — R 


( 2 . 22 ) 


DMj 


for r sa i?DM- Substituting this approximation for V into eq. (2.16c) and using 


eq. (2.20), we find that for r ss Rdm, F{r) approximately solves the differential 
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F(r) 


• 10-7 



Figure 2.4: Five scalings of the n = 1 excited state from fig. 2.8 (shown in blue) 
using a = 1.1,1.2,1.3,1.4,1.5 in the scaling formulas (2.21). Here we only show the 
dark matter prohle function F{r). 


equation 


^ , 2^ 2T2M(i?DM),^ „ 

Frr F Fj. „ (r F. 

r Rhu 


Multiplying through by r, we can write this differential equation as 

2T2M(i?DM), 


{rF)rr = 


-^DM 


-(r - RDM)irF). 


(2.23) 


Compare this to the Airy differential equation y" = xy whose general solution is 
a linear combination of the two standard functions Ai(x) and Bi(a:). Both Ai(x) 
and Bi(a:) oscillate for negative x and have exponential behavior for positive x. The 
function Ai(a:) is shown in £g. 2.5 It exponentially decreases to zero for a: > 0, 


whereas the function Bi(x) (not shown) exponentially increases for a: > 0. 



Figure 2.5: In black, the standard solution Ai(a:) to the Airy differential equation 
y" = xy. In brown and dashed, the approximation (2.28) to Ai(x) for x < 0. 
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Looking at eq. (2.23), we see that for r % Rbm, approximately solves an 
equation of the form y" = a{x — xo)y. The general solution to this differential 
equation is A • Ai(a^/^(a; — xq)) + i?-Bi(a^/^(a; —xq)). Since F exponentially decreases 
for r > Rbm, we must have, therefore. 


F(r) a - Ai 


2T^M{Rbm)V^'' 


-^DM 


-Rdm) 


(2.24) 


for r sa Rbm- If we choose the proportionality constant in eq. (2.24) so that the 
approximation has the correct value at r = i?DM, then we get 


F(r) % F{Rbm) 


R 


DM 


r Ai(0) 


Ai 


2T^M{Rbm) 


-^DM 


1/3 


(r — Rbm) 


(2.25) 


for r % i?DM- 

This expression contains the four numbers T, Rbm, M{Rbm), and F(i?DM)- 
These cannot all be independently chosen, for the static states effectively live in 
a 3-parameter space. Thus, we should be able to hnd some relation between these 
four numbers. In fact, we will argue that 


\F{Rbm)\ ^ for large n (2.26) 


where 


2-17/12^1(0) 0.133 


(2.27) 


Our goal is to derive eq. (2.26) beginning with the approximation (2.25). We will 
utilize a well-known approximation for the Airy function Ai(a:) for negative inputs: 

^2^3/2 + I) 


, , sin (fa:'" 

Ai(-x) ^ 


TTX 


1/4 


for a: > 0. 


(2.28) 


This approximation can be seen in £g. 2A It is extremely good for x > 1. We will 
also use the following dehnite integral, which is easily obtained: 


sm2 (|x3/2-f f) 

^-^dx = + 0(1). 


X 


1/2 


(2.29) 


Using eqs. (|2.16a|) and (|2.25|), we have 
M{Rbm) = I dr 


rRoM 

Stt 

Jo 


F{R 


DM 


)2i?2 


DM 


0 Ai(0)2 

F{Rbm)^RIm 


Stt- 


Ai(0)" 


Ai 


Ai 


2T^M{Rbm] 

-^DM 

2T^M{Rbm] 

p2 

r^DM 


1/3 


n 2 


(-Rdm — r ) 


dr 


1/3 


H 2 


dr. 
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Making the subsitution x = 
obtain 


We will temporarily abbreviate -Rdm = R, = M, = F, Ai(0) = A. 

2 \ 1/3 _ 

2 T M \ ^ using the approximation (2.28), we 


Ai(—a;)^ dx 

2 /2^3/2 ^ 


~A^ 


M ss Stt 


2T2]^^\“h3 r{2T^MRy/^ 


R^ 


98/3 p[ 2 T 2 MK]V 3 

4-T-A3F2i?A3jp^-i/3 


sm“ (|a; 


X 


1/2 


dx. 


Using the dehnite integral (2.29), we obtain 

28/3 


M 


A2 

217/6 


A2 


T-2/3^2^8/3j\^- 1/3 [(2T2Mi?)^/® + 0(1)] 

'Y'^l/3^2^17/6j^^-l/6 


(We dropped the 0(1) term because the other term dominates for large n.) Rear¬ 
ranging, we get 


Taking square roots gives us eq. (2.26), which completes the argument. 



0 100 200 300 400 500 600 700 800 


Figure 2.6: The ratio given in eq. (2.30) plotted as a function of n (n = 
0,1,..., 800). This graph shows that the approximation (2.26) is very good for 
large n and not too bad for small n. 


One might wonder how good the approximation (2.26) is for small n. In £g. 2.6 
we have graphed the ratio 


_|F(Rdm)|_ 

2-17/12 Ai(0) . TV6 M(Rdm)"/'2/?^W12 


(2.30) 


as a function of n. For the ground state (n = 0), the ratio is approximately 1.15, 
which is reasonably close to 1, and for large n the ratio is extremely close to 1. 
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Combining eqs. (2.25) and (2.26) gives the following approximation for the F{r) 
function of a static state of order n. (By convention F(0) is positive so the factor 
(—1)’^ is necessary.) 


-l)”2-l"/12T2/12M(/2DM)'/l'i?DM ' • J Ai 

Y2T2M(/2dm) 

.1/3 

j (r — T^dm) 

V ^DM 



(2.31 


This approximation is shown in hgs. 2.7 to 2.11 
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Figure 2.7: A ground state {n = 0) with T = 101y“^. Units are geometrized and 
given in light-years. The hrst three functions graphed are the solutions M(r), V{r), 
F{r) to eq. (2.16). They are, respectively, the total mass prohle, the Newtonian 
potential, and the dark matter profile. The fourth function is a rotation curve v{r) 
calculated in accordance with eq. (3.2). The vertical gray line marks the location 
of Rbm 


The dashed curves show the approximations described in section 2.5 


m 


eqs. (2.22) and (2.31). 
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Figure 2.8: An excited state (n 
description. 


• 10 ^ 


1). See the caption to fig. 2.7 for a fuller 
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An excited state (n = 2). See the caption to fig. 2.7 for a fuller 
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Figure 2.10: An excited state (n 
description. 


• 10 ^ 


10). See the caption to fig. 2.7 for a fuller 
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• 10 ^ 


Figure 2.11; An excited state (n = 100). See the caption to fig. 2.7 for a fuller 
description. The top and bottom of the F{r) function have been cropped out to 
show more detail. 
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3 


Wave Dark Matter and The Baryonic Tully-Fisher 

Relation 


In this chapter we describe some possible applications of the mathematical results 
from the previous chapter to theoretical astrophysics. 


3.1 Introduction 


In section 1.3 we introduced the theory of wave dark matter whose underlying equa¬ 


tions were the Einstein-Klein-Gordon equations (2.1). In this chapter for definiteness 
we fix a value for the fundamental constant T. We use T = 101y“^, 
tent with constraints from other studies 


17, 7, 45 


which is consis- 
of wave dark matter. 

We are interested in studying dark matter in disk galaxies and, for reasons out- 
possible connections to the baryonic Tully-Fisher relation. 


lined in section 1.4.4 


Although in disk galaxies the regular matter mostly lies in a disk, the dark matter is 
thought to lie in a “halo” which is at least approximately spherically symmetric. Sin 
proposed in that each galactic dark matter halo corresponds (in a first approxi¬ 
mation) to one of the spherically symmetric static states described in chapter If 
we take this suggestion as a starting point, a natural concomitant idea is that the 
baryonic Tully-Fisher relation might arise because a corresponding “Tully-Fisher- 
like” relation holds for the wave dark matter halos. This Tully-Fisher-like relation 
would hold because of the nature of wave dark matter. The fundamental question 
we ask in this chapter is whether the equations of wave dark matter have the po¬ 
tential to give rise to such a Tully-Fisher-like relation (see Question]^. We answer 
this question in the affirmative and exhibit two boundary conditions, which could 
be physical conditions imposed at the edge of each wave dark matter halo, that give 
rise to just such a Tully-Fisher-like relation. The boundary conditions are as follows: 


BCl: Fixing a length scale at the outer edge of halos implies a Tully-Fisher-like 


33 









relation with slope 4. 

BC2: Fixing the dark matter density at the outer edge of halos implies a Tully- 
Fisher-like relation with slope 3.4. 

In the rest of this chapter we work with the static states of the Poisson-Schrodinger 
system as models for dark halos. We can then freely use the scalings described in 
theoremprovided we are careful to stay in the low-field, non-relativistic limit. Our 
results will then be valid for the Einstein-Klein-Gordon system as well. 


3.2 Scaling Conditions 

We remind the reader of the discussion at the end of section [2~4| regarding the scalings 
of the static states. With T fixed, one more condition is sufficient to “fix a scaling” 
and give a unique sequence n = 0,1, 2,... of static states. We ask: 

Question 1. Are there are any scaling conditions for which the sequence of static 
states n = 0,1,2,... obeys a Tully-Fisher-like relation 

M , , 

— = constant (3.1) 

yX 

for some 3 ^ x ^ 4? 


We need to make this question more precise, as it is unclear above what M and v 
are, exactly, given a particular static state. Astrophysicists attempting to calibrate 
the baryonic Tully-Fisher relation are faced with a similar problem. In that case, M 
is an estimate for the total baryonic mass. For n, some take it to be the maximum 
observed velocity; some take it to be an average velocity computed in some well- 
dehned way; some take it to be the velocity observed at a particular (arbitrarily 
chosen) radius. It is well-known that the rotation curves of many spiral galaxies are 
somewhat flat, and in this case any of these three choices will produce similar values 
for V. 

We have chosen to take M = and v = v^Rdm)- This latter quantity is 

the circular velocity which would be observed for objects orbiting at r = Rbm, and 
is computed in the standard Newtonian way: 


accel. = 


n(r)^ M(r) 


v{r) = 


M{r) 


(3.2) 


Perhaps a more natural choice for M would be 


00 ? 


i.e. 


the total mass of the static 


state. But M(i?DM) Moo since i?DM is by dehnition the “outer edge” of the halo 
(see dehnitionand hg. |2.2[ ) , so it does not really matter which we choose. Similarly, 


another natural choice for v would be Umax, fhe maximum value of the function v{r) 


But again, by examining hgs. 2.7 to 2.11, one can see that u(i?DM) Umax- Thus, 


the following more precise version of Question [T] seems reasonable: 
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Question 2. Are there are any scaling conditions for which the seqnence of static 
states n = 0,2,... obeys a Tnlly-Fisher-like relation 


M{R 


dmJ 


v{R 


= constant 


dmJ 


(3.3) 


for some 3 < a; ^ 4? 


Now, in one sense the answer to Qnestion is a trivial “yes”, for we can jnst 
pick any x we like and let eq. (3.3) itself be the scaling condition. However, this is 


obvionsly cheating; we are interested in scaling conditions which might correspond 
to something physical. We will consider a scaling condition to be physical if it hxes 
some property of the dark matter prohle fnnction F{r) either at a particnlar point (a 
local condition) or in the large (a global condition). Orthogonally to the local/global 
distinction, we will consider scaling conditions which £x a horizontal property of F 
or a vertical property of F. The hrst type of scaling condition hxes some type of 
length, whereas the second type of scaling condition fixes |F| in some way. Looking 
at eq. (2.16a), we see that 2F^ gives mass density so that hxing a vertical property 


of F corresponds to hxing the density of the dark matter in some way. 

For local conditions, there are really only two natnral places to impose a scaling: 
at r = 0 or r = i?DM- We consider both, with horizontal and vertical scalings. In 
the category of global conditions, a horizontal scaling and a vertical scaling give two 
more scaling conditions for a total of six scaling conditions in all. We will discnss all 
these conditions in sections 13.41 to 13.61 after we make the dehnitions in section 13.31 

3.3 Definitions 

Definition 3. The F{r) function for any static state has an exponentially decreasing 
“tail” lying to the right of Rbm- There is some R > Rbm such that F{R) = ^F{Rdm); 
we refer to R — Rbm as the halflength of the state. 

Definition 4. For a Poisson-Schrodinger static state, the wavelength of F{r) at 
0 < r < Rom is given by the expression 

271 


2T{T -u + TV{r))' 


This definition follows naturally from an examination of eq. (2.16c). 


Definition 5. In the solution f{t,r) = F(r)e*“^* (see eq. (2.5)), the constant u is 


like a frequency but it is not a physical quantity because t is just a coordinate. We 
define 

UJtrneir) = (3.4) 

This is a physical quantity—the true frequency of the dark matter that would be 
measured at a particular value of r. The factor comes from the metric (2.4). 
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Definition 6. The average wavelength of a static state of order n is defined to be 


R 


DM 


n/2 + 1/4 


The reader is invited to examine figs. 2.7 to 2.11 to see why this is a sensible definition. 


Definition 7. The average density of a Poisson-Schrodinger static state is defined 
to be 

^ dvrr^ . 2F(r)2 dr. (3.5) 

Jo 


-ttR^ 

3 " -^DM 




DM Jo 


3.4 Boundary Conditions at r = B 


DM 


In this section we discuss imposing local scaling conditions at r = i?DM, both hori¬ 
zontal and vertical. For obvious reasons we refer to these two scaling conditions as 


boundary conditions, BCl and BC2. These were given in section |3.1| along with 
the Tully-Fisher-like relations they imply. We restate them here more precisely: 

BCl: Fixing a length scale at Rdm implies a Tully-Fisher-like relation with slope 4. 

BC2: Fixing |F"(i?DM)| implies a Tully-Fisher-like relation with slope 3.4. 

Note that, as we stated above, fixing |F(i?DM)| is equivalent to fixing the dark matter 
density at i?DM- 


d.^.i Boundary Condition 1: Fixing a Length Scale 

With regards to BCl, we should explain what “fixing a length scale” means. One 
easy way to understand this concept is in terms of the halfiength of static states—see 
definition Requiring the halfiength of all states to be the same is one way to fix a 
length scale at Rdm- Alternatively, we can fix the true oscillation frequency of the 
dark matter, a;true(’") (see definition]^, at a fixed distance outside (or inside) Rdm- 
This condition is explored more extensively in 1^. But most precisely, the best way 


to understand what fixing a length scale means is via the approximate solution ( |2.24[ ) 
from chapter which we repeat here. We discovered that for r 


Rdm, 


F(r) a - Ai 
r 


2T^M(Rdm) 

-^dm 


1/3 


(r — Rdm) 


(3.6) 


Up to a vertical scaling, the shape of R(r) around r = Rdm is given by the 


expression on the right side of (3.6). Fixing a length scale around r = Rdm evidently 


amounts to fixing the value of the constant 

2T^M(Rdm) 

-^DM 


(3.7) 
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This explains why fixing a length scale is basically eqnivalent to fixing a valne for 
the halfiength or to any of a nnmber of other conditions (snch as the one involving 
i^tme(’") stated above) which fix some length aronnd -Rdm- All snch conditions are 


fixing a particnlar horizontal scaling of the expression in eq. (3.6) 


Now, the rotation curve calculation (|3.2) shows that 




DM 


^{Rdu)^ ^ ^ 

Since T is a constant, we see that fixing a length scale is equivalent to fixing the 


constant in eq. (3.8). This explains why the scaling condition BCl implies a Tully- 


Fisher-like relation for the static states with slope x = A. 

3 . 4.2 Boundary Condition 2: Fixing the Density 

The truth of BC2 follows from eq. ( |2.26 ) from chapter]^ which we repeat here: 

|F(/?dm)| 2-iAi2Ai(0) . for large n. (3.9) 


From eq. (3.9) we see that fixing |F(i?DM)| amounts to fixing M^RdmYRdm- 


eq. (^), Rdm = M{Rdm)v{Rdm) so fixing M(i?DM)^.RDM is equivalent to fixing 

Af(/2DM)“^°'y(i?DM)^'^ = (Af(i?DM)/'y(-RDM)^'"^)“^°, i-o. equivalent to fixing 

M{Rdm) 

v{Rdm)^'^ 

Thus BC2 implies a Tully-Fisher-like relation with slope x = 3.4. 

3 . 4.3 Numerical Evidence 

In addition to the theoretical arguments we have just adduced for the truth of BCl 
and BC2, we can also check numerically that BCl and BC2 do give Tully-Fisher- 
like relations with the stated slopes. We used a computer to derive particular static 
state solutions for n = 0,1,..., 800 and then scaled these solutions in accordance 


with eq. (2.21) so that, first, they all had the same halfiength, and second, so they 
all had the same value for |F(i?DM)|- (The precise values for the halfiength and 
|.F"(.Rdm)| were chosen so that the resulting static states would be of reasonable sizes 
to be representing galactic halos.) The results are shown in fig. 3.1, in which we have 


plotted M{Rdm) vs. n(i?DM) in log-log space. For large n, the static states form a 
line with the stated slope. (For small n, the states deviate slightly from the correct 


slope because the approximation eq. (2.25) is not as good for those n. 


3.5 Central Conditions at r = 0 

In this section we discuss imposing local scaling conditions at r = 0, both horizontal 
and vertical. For obvious reasons we refer to these two scaling conditions as central 
conditions, CCl and CC2. 


37 
















BCl: Fixed Halflength 


BC2: Fixed |F(/2 dm)| 



^’(-Rdm) (km/s) ^(-Rdm) (km/s) 

Figure 3.1: Numerical evidence for BCl and BC2 (see page 36). At left, the static 
states n = 0,1,... ,800 with T = 101y“^ and a fixed exponential decay halflength 
of 884 ly. The dashed line has slope 4. At right, the static states n = 0,1,..., 800 
with T = 101y“^ and a fixed value for |F(i?DM)| of 2.473 x 10“®. The dashed line 
has slope 3.4. In these plots we have chosen not to adhere to our convention of using 


geometric units to facilitate comparison to fig. 1.3 


CCl: Fixing the wavelength of F at r = 0 (variant: fixing the location of the first 
node for n ^ 1) does not give a Tully-Fisher-like relation. 

CC2: Fixing F(0), i.e. fixing the density of the dark matter at the center of halos, 
does not give a Tully-Fisher-like relation. 


We do not provide theoretical arguments here for the truth of CCl and CC2, only 


numerical evidence. See fig. 3.2 


3.6 Global Conditions 

In this section we discuss imposing the following two global scaling conditions, the 
first a horizontal condition, the second a vertical condition. 

GCl: Fixing the average wavelength of F does not give a Tully-Fisher-like relation. 
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CC2: Fixed F(0) 


CCl: Fixed Central Wavelength 
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^ 10^° 
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^^(-RDM) (km/s) v{Rbm) (km/s) 

Figure 3.2: Numerical demonstration that CCl and CC2 do not give Tully-Fisher- 
like relations. At left, the static states n = 0,1,..., 800 with T = 10 ly“^ and a hxed 
central wavelength of 1500 ly. (Alternatively, hxing the location of the hrst node for 
n ^ 1 gives a qualitatively similar plot.) At right, the static states n = 0,1,..., 800 
with T = 101y“^ and a hxed value for F(0) of F{0) = 8 x 10“^. The dashed lines 
have slope 4. In these plots we have chosen not to adhere to our convention of using 


geometric units to facilitate comparison to hg. 1.3 


GC2: Fixing the average density of halos gives a Tully-Fisher-like relation with 
slope 3. This result is purely a Newtonian gravitational result and does not 
depend on the nature of wave dark matter. 


We do not provide a theoretical argument here for the truth of GCl, only numerical 
evidence. See hg. 3.3 The truth of GC2 follows from a simple Newtonian argument 
which does not use the specihcs of the wave dark matter model or the static states. 
Fixing the average density (see dehnition]^ amounts to hxing M(Rj^m)/R^m- The 
standard rotation curve calculation ( 3.2[ ) gives Rbm = M(Rbm)/v{Rdm)‘^■ Substitut¬ 
ing in, we see that hxing the average density amounts to hxing n(i?DM)®/Af(i?DM)^, 
or equivalently, hxing 


M{Rbm) 

v{Rdm)^ 


Thus we get a Tully-Fisher-like relation with slope x = 3. The argument just given 
is basically the same as the one given to support a Tully-Fisher relation with slope 3 
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in the ACDM paradigm (see 50, O), with Rum replaced by, for example, r 2 oo, the 


radius enlosing a region with average density equal to 200 times the critical density 
of the universe. 


GCl: Fixed Average Wavelength 
10^3 


GC2: Fixed Average Density 




(km/s) 

Figure 3.3: Numerical demonstration that GCl does not give a Tully-Fisher-like 
relation and GC2 does. At left, the static states n = 0,1,..., 800 with T = 10 ly”^ 
and a hxed average wavelength of 1500 ly. At right, the static states n = 0,1,..., 800 
with T = 10 ly“^ and a hxed average density of 10“^^. The dashed lines have slope 
3. In these plots we have chosen not to adhere to our convention of using geometric 


units to facilitate comparison to hg. 1.3 


3.7 Discussion and Conjectures 

In the previous section we discovered that to get a TuIIy-Fisher-like relation in the 
context of wave dark matter we need to impose a boundary condition (BCl or BC2). 
We could also impose the global condition GC2, but this does not use any of the 
characteristics of wave dark matter. Because GC2 applies to every theory of dark 
matter, we will ignore it in the discussion which follows. 

3.7.1 The Slope of the Baryonic Tully-Fisher Relation 

In this section we comment briehy on the connection between a Tully-Fisher-like 
relation for wave dark matter halos and the observed baryonic Tully-Fisher relation. 
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Define the haryon fraction of a galaxy to be the fraction of the galaxy’s mass which 
is baryonic. If the baryon fraction were the same for every galaxy, then the baryonic 
Tully-Fisher relation would have the same slope as a Tully-Fisher-like relation for 
dark halos. However, it is observed that the baryon fraction of galaxies decreases 
with size. Small galaxies contain mostly dark matter, and only very large galaxies 
contain baryons in sufficient quantities to approach the cosmic baryon fraction, which 
is around 1/6. This is known as the missing baryon problem (see ^f\). 

In M vs. V log-log space, when switching from plotting dark mass to baryonic 
mass, all galaxies shift vertically downward. The missing baryon problem implies 
that large galaxies shift less than small galaxies. Therefore, if the baryonic Tully- 
Fisher relation is related to a Tully-Fisher-like relation for dark halos, the former 
relation should have a steeper slope than the latter. 

3. 7 .2 Conjectures 

In this section we describe a testable prediction of the wave dark matter model. 
Assuming a new idea called “dark matter saturation” described below and made 
precise in Physics Conjecture given the distribution of the baryonic matter in 
a galaxy, we should be able to compute the distribution of the wave dark matter 
and hence the total mass distribution, as long as everything is approximately static 
and spherically symmetric. With the total mass distribution in hand it will be 
possible to compare to observations—for example, to compute rotation curves. These 
predictions will be possible once specific stability questions, described below, are 
answered. 

If we assume that dark matter in many galaxies is approximately static and 
spherically symmetric, then it makes sense to look for static, spherically symmetric 
solutions to wave dark matter. Referring back to the scaling equations, we see that 
there is a two-parameter family of static spherically symmetric solutions. These two 
parameters are (n, a), where n ^ 0 is an integer referring to the excited state and 
a is the scaling factor. Recall that n = 0 refers to the ground state, and that more 
generally n is the number of zeros of F{r). 

Now suppose that we solve for wave dark matter solutions in the presence of 
regular matter (which, for our purposes here, we will also assume is spherically 
symmetric and static). While the regular matter, through gravity, will change each 
wave dark matter solution, we still expect to find a two parameter (n, a) family of 
solutions. 

One great beneht of our discussion so far is that it removes one of the parameters, 
namely the continuous parameter a. In this chapter, we showed that a boundary 
condition (BCl or BC2) is roughly what is needed to recover a Tully-Fisher-like 
relation. To be clear, we have not explained why such a boundary condition should 
be expected from the theory, just that something close to it seems necessary to be 
compatible with the observations which make up the baryonic Tully-Fisher relation. 
In any case, assuming one of these boundary conditions effectively determines a, 
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leaving only one parameter free, namely n. 

What should the value of n be? Given a regular matter distribution, for each n 
we get a precise wave dark matter distribution which satishes the above boundary 
condition. Some of these wave dark matter solutions will be stable and some will 
be unstable. Numerical results show that ground states of wave dark matter are 
stable while excited states, without any other matter around, are unstable 22 . On 
the other hand, 22 also shows that excited states may be stabilized by the presence 


of another matter field. Our conjecture is that the regular, visible, baryonic matter 
stabilizes wave dark matter in galaxies. In fact, given a regular matter distribu¬ 
tion, we conjecture that there exists a largest value of n, call it iV, for which the 
corresponding wave dark matter solution is stable. We conjecture that galaxies are 
described best by choosing n = N. 

The total mass of the spherically symmetric static states described in this paper 
increases with n and we expect the same to be true for the distributions of dark 
matter we are describing now. Thus, setting n = N is consistent with the idea that 
galaxies are “dark matter saturated”, meaning that they are holding as much dark 
matter as possible, subject to the boundary condition above. Since galaxies typically 
exist in clusters which are mostly made of dark matter, it seems likely that they are 
regularly bombarded by dark matter, so that it would be natural for them to reach 
this state of saturation n = N. 

To make this discussion precise, we need a model for the regular matter. In order 
to study stability questions, we need to know how the regular matter distribution 
changes as the wave dark matter distribution changes, and vice versa. 

For example, a relatively simple way to model regular matter is with another 
scalar held. There are others ways to model regular matter which we do not discuss 
here. We caution the reader that this second scalar held is only a practical device 
for approximately modeling the regular baryonic matter—namely the gas, dust, and 
stars in a galaxy. In no way are we suggesting a second scalar held should exist 
physically. Furthermore, the parameters of this second scalar held are chosen simply 
to ht the regular matter distribution of a galaxy as well as possible. 

Let /i exactly model wave dark matter with its fundamental constant of nature 
Ti. Let /2 be a convenient device for approximately modeling the regular baryonic 
matter consisting of the gas, dust, and stars of a galaxy, where T 2 , which is not a 
fundamental constant of nature, is chosen as desired to best ht the regular matter. 
The action is then 


Hg JiJ2) = 


Rr,-2A- 1671 




+ 1/2!' + 


\df2 


'Y'2 

2 


dVg, (3.10) 


where A is the cosmological constant and may as well be assumed to be zero for 
our discussion on the scale of galaxies. The above action results in the following 
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Euler-Lagrange equations: 


G + Ag = 8tt 


dfi (X) dfi + dfi (X) dfi 
T2 


i/ir + 


2 , M/ll 


T2 


+ 


df2 (8) ci/2 + df2 (X) df2 

Ti 


1/2/ + 


2 , M/2I 


'Y'2 

2 


□/i = T^/i 

□ /2 = Tlf2 


(3.11a) 


(3.11b) 

(3.11c) 


We approximate the regular matter distribution with a ground state solution for / 2 . 
We have two free parameters with which to approximate the given regular matter 
distribution, namely T 2 and the “scaling parameter” for the ground state solution, 
which we could call 02 . This should allow us to choose two physical characteristics of 
the regular matter. We choose to specify the total mass Mb and the radius Rb of the 
regular matter, perhaps defined as that radius within which some fixed percentage 
of the regular matter is contained. 

As described already, we impose the boundary condition BCl or BC2 for /i. 


Then for each choice of n ^ 0, we get a solution to the system of equations (3.11) 
which reduces to a system of ODEs in a manner very similar as before. Some solutions 
will be stable and some will be unstable. Since we now have the dynamical equations 


(3.11), these stability questions are now fairly well dehned. 


Hence, for each Mb, Rb, and n, we get a static, spherically symmetric solution to 
eq. (3.11) satisfying the boundary condition BCl or BC2. 


Math Conjecture 1. In the low-field, nonrelativistic limit, for each choice of total 
regular mass Mb and regular matter radius Rb, there exists an integer ^ 0 such 
that static, spherically symmetric solutions to eq. (3.11) satisfying the boundary 


condition BCl or BC2 with n ^ N are stable and those with n > N are unstable. 


If this math conjecture is true, or even if there is just a largest or most massive stable 
n, then there is a natural physics conjecture to make as well. 


Physics Conjecture 1 (“Dark Matter Saturation”). The dark matter and total 
matter distributions of most galaxies which are approximately static and spherically 
symmetric are approximately described by static, spherically symmetric solutions to 
eq. (3.11) satisfying the boundary condition BCl or BC2 with n = N. 


This last conjecture only leaves two parameters open, namely T, the fundamental 
constant of nature in the wave dark matter theory, and the length scale (or halflength) 
from BCl or the boundary dark matter density from BC2. Hence, there are effec¬ 
tively only two parameters left open with which to £t the dark matter and total 
matter distributions of most of the galaxies in the universe. Therefore the physics 
conjecture stated above should be a good test of the wave dark matter theory. 
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Appendix A 

The Einstein-Klein-Gordon-Maxwell Equations 


In this appendix we derive the Einstein-Klein-Gordon-Maxwell equations from an 
action using a variational principle. On a spacetime with metric g, complex scalar 
field /, and one-form A, we define the action 





2A+Iff+^-\dA\^ 


dV. 


(A.l) 


Here U is any precompact open set with smooth boundary, so that the integral 
converges. We will require that g, f and A are critical points of the functional J^u 
for all possible choices of U and all variations compactly supported within U. In the 
three sections below we vary each of g, f, A independently to derive the Einstein 
equation, the Klein-Gordon equation, and Maxwell’s equations. 

We will make use of the Hodge star operator ★ and the codifferential operator S 
which is adjoint to d with respect to the scalar product (•,•). For brevity we will 
omit the domain of integration U. Because the variations we consider are compactly 
supported, all boundary terms are zero. 


A.l Varying the Metric g 

In the calculation below we will use the following well-known results concerning a 
variation g = h of the metric: 

R = —(h, Ric) -f- div(div h) — □(tr h) 

dV = - (tr h) dV. 

We will also need to know the variations of \df'^ and \dA'f. First, since gx^g^’" = 
differentiation gives gx^g^^ = —hx^g^^- Multiplying by g^'^ gives g'^P = —R'p. To 
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calculate the variation of \dff, begin with 


W\" - g^“fxU. 


Differentiating, we get 

\dff = -h^^fxl, = -{K df ® df) = -<h, df^ df). 

Therefore 

\df\^ = -{K]^{df®dJ+dJ®df)). 

To calculate the variation of write F = dA so that in coordinates, 

\dA\^ = 


Differentiating, we get 

\dA\^ = -h^^^g'^PFx.F.p - g^^F'PFx.F^p. 

The two terms above are equal by the anti-symmetry of F. Thus 

\dA\^ = -2(h,Ci3{dA(g)dA)) 


where Ci 3 {dA(S)dA) denotes the tensor obtained by contracting on the hrst and third 
slots of dA (X) dA. 


With these results in hand, if we perform the variation in eq. (A.l), we get 
0 = 


o- \ , o / r, df ®df + df ®df ^ j 

- <h, Ric> -t- 87r( h, -—-h Ci 3 {dA (x) dA) 


+ (tr h) 


- - A - Stt 
2 


|2 




dV. 


We immediately dropped the terms div(divh) and □(trh) because using the diver¬ 
gence theorem, they give boundary terms which are zero. If we note that tr h = (h, g'), 
the equation above takes the form 


0 


^-{h,Z)dV 


where Z is the (0, 2)-tensor 


Z 


Ric ~ 2 ^g F Ag — Stt 


df ^ df + df ^ df 


-I- Ci3{dA (X) dA) 





9 
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This equation must hold for all variations h, which implies that Z = 0. Thus the 
variational principle has led us to the Einstein equation 


G + Ag = SttT 


(A.2) 


(recall the dehnition G = Kic —{l/2)Rg) where the energy-momentum tensor T is 
given by 


T 


df ^ df + df ^ df 


'^2 



g + Ci;^(^dA (X) dA^ 



(A.3) 


A.2 Varying the Scalar Function / 

It suffices to consider the variation of the integral 

J [M/I" + T"|/|"] dV- ( [(df,df) + TV/1 /I"' 

Performing the variation, we get 

0 = J [(df, df) + (df, df) + TV/ + TV/] dV 

= J[/(<5d + T2)/ + /(5d + T2)/] dV. 

This equation must hold for all variations /. Since we can vary the real and imaginary 
parts of / independently, this implies that {6d + T^)/ = 0 and {6d + T^)/ = 0. The 
second equation is the complex conjugate of the hrst and is superfluous. The operator 
Sd acting on functions is equal to 6d -I- dS, which is the Laplace-de Rham operator, 
the negative of the Laplace-Beltrami operator (in this case the d’Alembertian) □. 
Therefore the variational principle has led us to the Klein-Gordon equation 

□/ = TV- (A.4) 


A.3 Varying the One-Form A 

It suffices to consider the variation of the integral 

^\dAf dV = ^ dA A *dA. 

Performing the variation, we get 

0 = r dA A *dA + dA a *dA. (A. 5) 
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The two summands above are equal because, by the definition of the Hodge star, 

dA A *dA = dA A *dA = {dA, dA)dV. 


Since 


eq. (A.5) becomes 


d{A A -kdA) = dA A -kdA — A a d*dA 


0 = A A d*dA. 


This equation must hold for all variations A, which implies that d*dA = 0, and this 
equation is equivalent to -kd-kdA = 0. Recalling that, up to a sign, ★d* is equal to 
the codifferential operator h, the equation we derive from varying A is 


5dA = 0. 


(A.6) 


This is half of what are usually referred to as Maxwell’s equations. A more conven¬ 
tional way of writing them is to define the Faraday tensor F = dA. We automatically 


have dF = 0; this equation paired with eq. (A.6) becomes 


dF = 0 
5F = 0. 


A.4 The Einstein-Klein-Gordon-Maxwell Equations 

Requiring the metric the complex scalar function /, and the one-form A to be crit¬ 


ical points of the functional (A.l) has led us to the Einstein-Klein-Gordon-Maxwell 


equations collected in eqs. (A.2) to (A.4) and (A.6), which we repeat here: 


G + Ag = SttT 
□/ = TV 
5dA = 0. 


(A. 7a) 
(A.7b) 
(A.7c) 


where 


T = 


df ® df + df ® df 
'^2 


^ + \ f\j 9 + Ci 3 (dA (X) dA) — -\dAf g. 
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Appendix B 


The Einstein-Klein-Gordon Equations in Spherical 

Symmetry 


In this appendix we derive the system of ordinary differential equations for the spheri¬ 
cally symmetric static states of wave dark matter. The calculations we carry out here 
are a special case of calculations carried out in 
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we assume from the beginning that the metric (B.2) is static. 


but are much simpler because 


The spherically symmetric static states are solutions to the Einstein-Klein-Gordon 
equations 


G = 8tt 

□/ = TV- 


df ®df + df ®df (\dff , I ,,2 




■^2 


+ \fn9 


We begin with the spherically symmetric static metric 

g = — dt"^ + (l - (^ 7-2 ^2 ^q2 ^2 g|j^2 g ^ 


We define 


<h(r) = 1 


2M(r) 


so that the metric (B.2) can be written as 


(B.la) 

(B.lb) 

(B.2) 

(B.3) 


g =dt'^ + ^{r) ^ -1--I-sin^ (B.4) 
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Lemma 2. The nonzero Christoffel symbols of the metric (B.4) in the {t,r,6,(j)) 
coordinates are as follows: 


= -r$ 

= -r sin^ e $ 


r® = = r~^ 

^ rO ^ Or 

r® ^ = — sin 6 cos 9 

^ rcj) ^ (f>r ' 

rt - ri - cote. 


Proof. These results follow immediately from short computations using the usual 
formula for the Christoffel symbols: 


^A/i 2^ + 5 'pp,A fi'Ap,p)- 


□ 


Lemma 3. The nonzero components of the Ricci curvature tensor associated with 
the metric (B.4) in the (t,r,9,(j)) coordinates are as follows: 


Rictt = {Vrr + + 2r-^Vr)e^^^ + {l/2)Vre^^^r 

Ric,, = -{Vrr + 

Ricpe = 1 — $ — (l/2)r$r — rVr^ 

Ric00 = sin^ 6* (1 — <h — (l/2)r<hr — 


Proof. The formula for the components of the Ricci curvature tensor in terms of the 
metric and the Christoffel symbols is 


RicAp = - rL,p + .,, 


_ pP -pi' 

^ /ip* 


(B.6) 


In the calculations below we enclose each of the four terms in eq. (B.5) in square 


braces, omitting all the Christoffel symbols which are zero. Then we make whatever 
cancellations are possible, substitute in for the remaining symbols, and simplify. We 
have 


Wq. - - [] + [r:,(rf, + r;, + C + rf*)] - [U.r;, + r;,ra 

— -I- F'" F^ -i- F'’ F® J- F'' F*^ _ F* F'" 

tt,r ' tP rr ' ^ r9 ' ^ rcj) ^ tr^ tt 

= {Vre^^^)r - (l/2)R,e2^<h, + 2r"Wre^^^ - 
= Vrre^^^ + 2V^e^^^ + - (l/2)V;e^^<h^ + 2r-V,e2^<h - 

= {Vrr + v^ + 2r"Vpe2^4) + (l/2)Ke^'^<h. 


and 


Ricr - [rL,] - [] + D - [] = 0 
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and 


Rico -[]-[] + []-[]-•) 

and 


Ric* - D -[] + []-[]- 0 

and 


Ric„ - [r;„] - [r'„, + + r"s_, + + [r;,(rt, + r;, + r", + rf^)] 

_ rpt pt I pr pr , p^ p6» , p</> p</< i 
[i “T i “T i -r L 

_ _-pt _ p0 _ p0 I pr pt I pr -p9 I pr p0 _ pt pt _ p0 p0 _ p</> p0 

'rt,r 'r0,r r(/),r rr-*- rt ' ^ rr^ r9 ' ^ rr^ r(j) ^ rt^ rt r9^ r9 ^ r<l}^ r<l) 

= -Vrr + 2r-2 + (-(l/2)$-i$^)K + (-(l/2)<l>-^$^)(2r-^) - - 2r-^ 

= -{Vrr + V^) - {l/2)Vr^~^^r “ 


and 


Ric,« = [rj,„] - [rt„ + r;,,, + C,, + j + [Crjj - 

— /'F^ — F^ ^ 

^ 9(p\^ r9 ^ r<p) 

= cot 6{r~^ — r~^) 

= 0 


Ric,* - [rf„] - [rt,,^ + T%,^. + r“,_^ + t%J + []-[]- o 


Ric« = [r;,,] - + [r;,(r‘, + r', + C + t%)\ - [Cr;, + tiA + r^rtl 

_ pr _ p</> I F’’ r* -I- r*” I F»- F</' _ F?- F^ _ F*/* F*/* 

~ ee,r ^ 60,0 ^ ^ ee^ rt ^ ^ ee^ rr ^ ^ ee^ r<j) ^ ee^ er ^ e<j>^ e<f> 

= (-r$)^ - (cot 0)0 + {-r^)Vr + (-r$)(-(l/2)$”^$,) 

+ (—r$)r“^ — (—r<l>)r“^ — cot^ 6 
= —$ — r$r + csc^ 6 — rVr^ + (l/2)r<l>r — cot^ 6 
= !-<!>- {l/2)r^r - rVr^ 


and 

Ric« - VUJ - VtJ + []-[] = 0 
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and 


- [] + [r;^(rj, + r;, + r"„ + r*j + r%r*^,] 

rp(/> pr I ■p(l> -pO I -pr ■p(l> , -pO p(/> i 

(pA (P(p~^ ^ (pe^ (P(p ^ (p<p^ (Pr ^ (p(p^ (pel 

_ pr I -pO I pr -pt I pr pr , pr -pO pr p0 p0 p0 

^ (p(p,r ^ (p(p,e ' (p(p^ rt ' ^ (p(p^ ~r L ^.q i i 

= (—r sin^ 6 $),, + (— sin 9 cos 6)0 + (—r sin^ 9^)Vr 

+ (—r sin^ 6* $)(—(1/2)$”^$^) + (—r sin^ 0 d))r“^ — (—r sin^ 6* $)r“^ 

— (— sin 9 cos 9) cot 9 

= — sin^ 9^ — r sin^ 9^^ — cos^ 9 + sin^ 9 — rVr sin^ 9 $ 

+ (l/2)r sin^ 9^r + cos^ 9 

= sin^ 6* (1 — $ — {l/ 2 )r^r — rVr^). 


These calcnlations snffice becanse the Ricci cnrvatnre tensor is symmetric. □ 


Lemma 4. The scalar curvature associated with the metric 
coordinates is 


(B.4) in the {t, r, 6*, 0) 


R = -2{Vrr + Vp + 2r-V,)4> - + 2r-^{l - $ - r<h,). 


Proof. The scalar cnrvatnre is the trace of the Ricci cnrvatnre tensor. Using lemma|^ 
and the fact that the metric (B.4) is diagonal, 


R = RiCtt Ric,, Ric^^ +g‘^^ Ric^^ 

= {-e~^^){{Vrr + Vp + 2r-^Vr)e^^^ + (l/2)Ke2^<h,] 

+ ^-{Vrr + Vp) - (l/2)K<h"^<h^ - 
+ r~^[l - $ - (l/2)r<hr - rVr^] 

+ r~^ sin“^ 9 [sin^ 0 (1 — $ — (l/2)r$r — rl4-<h)] 

= -{Vrr + Vp + 2r-V,)<h - (l/2)K<h, - {Vrr + Vp)^ - (l/2)K<h, - 
+ 2r“^ — 2r“^<h — r~^^r — 2r“^14-<h 

= -2(U,, + Vp + 2r-V^)<h - Vr^r + 2r-2(l - $ - r<h,). □ 


Lemma 5. The nonzero components of the Einstein curvature tensor G associated 
with the metric (B.4) in the (t,r,9,(f)) coordinates are as follows: 


Gtt = r — <h — r^r) 

G„ = - $ - 2 rVr^) 

Gee = AiVrr + Vp + r-V,)<h + (l/2)r-^<h, + (l/2)U,<h,] 

sin^ 0 [(V;, + 1/2 ^ r-V,)4> + (l/2)r"i$, + (1/2)^$,]. 
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Proof. The Einstein curvatnre tensor is defined by eq. (1.2). (Remember that in this 
dissertation, A = 0.) Using lemmas and we have 


Gtt — Rictt — {1/2) Rgtt 

= [{Vrr + + 2r-V,)e2^$ + {l/2)Vre^^^r\ 

- (l/2)[-2(U,, + + 2r-V^)$ - Vr^r + 2r-2(l - $ - r$^)](-e2^) 

= — $ — r$j,) 


and 


Grr = RiCrr—{1/2) Rgrr 

= [-{Vrr + K") - {l/2)Vr^-^^r " 

- (l/2)[-2(U,, + + 2r-V^)<h - + 2 r-^{l - <h - 

= -r-2$“^(l-<l)-2rU.<h) 


and 


Gee — Ricee —{l/2)Rg0e 

= [!-$- (l/2)r<l>,. - rVr^] 

- (l/2)[-2(U,, + 1/2 + 2r“V^)<h - Vr^r + 2r-2(l - $ - r^r)y 
= ryVrr + 1/' + r-y)^ + (l/2)r-^<l>, + (1/2)U,$,] 


and 


G Ric,^(j/j {1/2) Rg^i^fj^ 

= [sin2 0 (1 - $ - (l/2)r$^ - rVr^)] 

- (l/2)[-2(U,, + + 2 r-y)^ - Vr^r + 2r-2(l - $ - r$^)]r2 sin^ 0 

= r^sin^^ [(R.^ + Vr^ + r-^)^ + (l/2)r-^<h, + (l/2)U,<h,]. 


All other components are zero 
is diagonal. 


by lemma and the the fact that the metric eq. (B.4) 

□ 


Einstein’s equation is G = SttT. For wave dark matter the energy-momentum 


tensor T is (refer to eq. (B.l)) 

df®df + df®df l\dy 


T = 


'f2 


'f2 


+ i/n9. 


(B.6) 


To solve the Einstein equation we need to equate the components of the Einstein 
curvature tensor with the components of the energy-momentum tensor, e.g. Gu = 
STrTtt, Grr = 8nTrr, etc. To do this we need to write down the components of T. 
We are solving the Einstein-Klein-Gordon equations in spherical symmetry, so / is 
a function of t and r only. 
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Lemma 6. With f = f{t,r), the nonzero components of the energy-momentum 
tensor (B.6) in the (t,r,9,(j)) coordinates are as follows: 


Ttt- 

Ttr = Trt = 

Trr - 

Tee = 

T(p(j) = 

= T~^ifJr + ftfr) 

= T-^r\-TW + e-^^\ftf-mf) 

= T-V sin^ 6 {-ryf + e-^^\f,f - «F|M^). 

Proof. First, we have 



\df\" - + Kl/rl" 


Thus 


r^Tu = 

^ ftft + ftft - {-e-^^\ftf + + r^\ff){-e^n 

-W^\f\^ + \ft\^ + e^^mf 


and 


= T^Trt = 

'■ ftfr + ftfr 


and 



^ frfr + frfr - + 0\fr\" + T"|/|")«>-‘ 

:-T=4-‘|/|" + e-^"'4>-‘|/d + l/rt 


and 


r^Tee = 

--{-e-^^\ftf + mf + TVfy 


and 


'y^2rTi 



The other components of T are zero because / depends only on t and r and the 
metric (B.4) is diagonal. □ 
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Theorem 7. In the spacetime with the static metric (B.4) and energy-momentum 
tensor (B.6), where f = f{t,r) is complex-valued, Einstein’s equation (B.la) reduces 
to the following four PDFs: 

1 - $ 


= Swr^llfl + T ‘e 

0 = ftfr + ftfr 
2ri xi2 


-2„-2y 


|/i|" + T-2$|/,|2] 


1 _ $ _ 2rVr<!> = Snr^[\ff - T~^e-^^\ftf - r-^<^\frf] 


TV2[(K. + + r-V,)$ 

+ (l/2)r-'$, + (1/2)K$.] = 87rr\-T^\ff + e-^^\ftf - <!>\frf). 


(B.7) 

(B.8) 

(B.9) 

(B.IO) 


Proof. Using the results of lemmasandset G = SvrT and write down the resulting 
PDEs. Note that equations coming from equating the 99 and 00 components are 
identical. □ 


In the Einstein-Klein-Gordon system, the Einstein equation (B.la) is coupled to 
the Klein-Gordon equation (B.lb). The following lemma is the partner to theorem]^ 

Theorem 8. In the spacetime with the static metric (B.4), where f = f{t,r) is 

DE: 

(B.ll) 


complex-valued, the Klein-Gordon equation (B.lb) reduces to the following PDE: 
^-^^ftt + Vr^fr + (l/2)<h,/, + 2r”l<h/, + ^frr = 


— e 


Proof. We use the well-known coordinate expression for the d’Alembertian 
where denotes the absolute value of the determinant of the matrix representing 


g. For the metric (B.4), \g\ = V^sin^^. Thus 

□/ = (e-^4)^/V-2sin-^0)[(e^$~i/Vsin0 • -f sin0 • ^fr)r] 

= sin-^ 9) sin 9 fu + sin 9 /, 

-I- (l/ 2 )e^<l)“^'^^<l)rr^ sin 9 fr + sin 9 fr + sin 9 f^r] 

= -e-^^ftt + Vr^fr + {l/2)^rfr + 2r-^^U + ^frr. 

The result follows. □ 

We have shown that to solve the Einstein-Klein-Gordon system in the metric 


( |B.4[ ) with / = /(t,r), it suffices to solve eqs. ( |B.7[ ) to ( |B.11[ ). This system of 
equations is overdetermined. We have the following theorem: 


Theorem 9. In the spacetime with the static metric (B.4), where f = f{t,r) is 
complex-valued, to solve the Einstein-Klein-Gordon system (B.l) it suffices to solve 
eqs. (B.7), (B.9) and (B.ll), which we repeat here: 


1 _ $ _ = 8nr^[\f\^ + 

1 _ $ _ 2rW<h = Snr^ilff - T-^e-^^\ftf - r~^<^\frf] 


(B.7) 


—e 


-2V 


fu + + (1/2)$,/, -f 2r-i$/, -f $/,, = TV- 


dBlTl) 
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Proof. Assuming eqs. (B.7), (B.9) and (B.ll) hold, we show eq. (B.8) holds. By- 
adding and subtracting eqs. ( B.7| ) and (B.9), we get 

1 - $ - (l/2)r<l>^ - rVr^ = Svrr^l/I^ (B.12) 

-(l/2)r<h, + rVr^ = 87rr2T-2[e-2^|/i|2 + $|/,|2]. (B.13) 


Differentiating eqs. (B.12) and (B.13) with respect to t, we get 




0 — [e + ^\fr\ ]t — e {fttft + fttft) + ^{ftrfr + ftrfr)- 


-2V 


(B.14) 

(B.15) 


Using eqs. (B.ll), (B.14) and (B.15), we have 

+ ftfr)r = ^{ftrlr + ftrfr) + ^{ftfrr + ftfrr) 

= {fttft + fttft) + ^{ftfrr + ftfrr) 

= (~e ftt + ^frr)ft + ( —+ *h/rr)/t 

= r\fjt + fft) - (14-$ + (1/2)$, + 2r-'$)(/,7, + frft) 

= -(K$ + (1/2)$, + 2r-i$)(/,7 + 7,7). 

Thus 

(77 + ftfr)r = -{Vr + 2r-l + (1/2)$-!$,) (77 + fjt). 

Using the equation just obtained, we then have 

[r7^$'/777 + 77)]. = 2re^$'/777 + 77) + rV,e'^$'/777 + 77) 

+ (i/2)r7^$-'/2^,(77 + 77) 

- 7e'"$V7U, + 2r-i + (1/2 )$-i$,)(77 + 77) 

= 0 . 

Therefore for each fixed t the function r^e^$^/^(77 + 77) is a constant. Taking 
the limit as r ^ 0 shows this constant must be zero. Si nce is nonzero for 

all r > 0, we conclude that 77 + 77 = 0) which is eq. (B.8). 


Now assuming eqs. (B.7) to (B.9) and (B.ll) hold, we show eq. (B.IO) holds. 
Differentiating eq. (B.9) with respect to r, we obtain 


$, - 2U,$ - 2rU,,$ - 2rU,$, = 


2-2V\ 


T-^$|7l1 + 87r7[7/ + 7/ + 2T-V,e-2^|7|" 


167rr[|/| — T 

- t-7-2^(7.7 + Irft) - T-2$,|7I' - t-2$(7,7, + 7,7)]. 
Multiplying through by —(l/2)T^r gives 

T77U,,$ + r-V,$ + U,$, + (l/2)r"^$,] = 

87r7[-T2|/|^ + e-2^|7|^ + 4>|7I^] - 47r7[T77/+ If) + 2U,e-2^|7|^ 

- e-^^{ftrft + Uft) - <^r\frf " $(7.7 + 7.7)]- (B.16) 
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Multiplying eq. (B.13) by T^rl4- gives 


TV(V;2$ - (l/2)K<^>.) = ATTr^[2Vre-^^\ftf + 2V;$|/.|1 


Adding eqs. (B.16) and (B.17), we get 


(B.17) 


T\^[{Vrr + + r-V^)$ + (l/2)r-^d), + (l/2)V;d),] = 

S7ir‘^[~T‘^\ff + - ^frf] - 47rr^[T2(/^/ + /^/) - + ftrft) 

- - 2VMfr\^ - ^frrfr + Ufr) " 4r-'4)|/,|']. 


Comparing with eq. (B.IO), we see that it now suffices to show that the last expression 
in square brackets above is zero. We can rewrite this expression as 


(TV - (l/2)<hV, - v;4>/. - ^frr - 2r-^^fr)fr 

+ (TV - (l/2)$Vr - Wr - 4>/rr- " 2r-l<h/,)/, 

- (ftrft+ Uft). (B.18) 


Using eq. (B.ll), this expression becomes 


.-2V/ 


■e (fttfr + fttfr) — e (ftrft + ftrft) — (ftfr + ftfr)t. 


This last quantity is zero as desired by eq. (B.8). This completes the proof. 
The goal of this appendix is the following theorem: 


□ 


iujt 


Theorem 10. In the spacetime with the static metric (B.2), where f(t,r) = F(r)e 
with F real-valued and u real, to solve the Einstein-Klein-Gordon system (B.l) it 
suffices to solve the following three coupled ODEs: 


Mr = 47rr2 . ^ [(T^ + + ^Ff] 


^2 
2 $ 


2 -2V 

uj e 


) - <hFV 

Frr + -Fr + K-F, + ~Fr = 4>-^ (T^ - F. 


(B.19) 

(B.20) 

(B.21) 


Proof. These three equations are rearranged versions of eqs. (B.7), (B.9) and (B.ll) 
from theorem]^ with the particular form for / substituted in and with eq. (B.3) 
used in places so that the function M(r) appears. □ 
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